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ABSTRACT 

We present an inventory of the cosmic mean densities of energy associated with ah the known 
states of matter and radiation at the present epoch. The observational and theoretical bases 
for the inventory have become rich enough to allow estimates with observational support for the 
densities of energy in some 40 forms. The result is a global portrait of the effects of the physical 
processes of cosmic evolution. 



Subject headings: cosmology 

1. Introduction 

There is now a substantial observational basis 
for estimates of the cosmic mean densities of all 
the known and more significant forms of matter 
and energy in the present-day universe. The com- 
pilation of the energy inventory offers an overview 
of the integrated effects of the energy transfers in- 
volved in all the physical processes of cosmic evo- 
lution operating on scales ranging from the Hub- 
ble length to black holes and atomic nuclei. The 
compilation also offers a way to assess how well 
we understand the physics of cosmic evolution, by 
the degree of consistency among related entries. 
Very significant observational advances, particu- 
larly from large-scale surveys including the Two 
Degree Field Galaxy Redshift Survey (2dF: CoUess 
et al. 2001), the Sloan Digital Sky Survey (SDSS: 
York et al. 2000; Abazajian et al. 2003), the 
Two Micron All-Sky Survey (2MASS: Huchra et 
al. 2003), the HI Parkes All Sky Survey (HIPASS: 
Zwaan et al. 2003), and the Wilkinson Microwave 
Anisotropy Probe (WMAP: Bennett et al. 2GG3a), 
make it timely to compile what is known about 



the entire energy inventory. 

We present here our choices for the categories 
and estimates of the entries in the inventory. Many 
of the arguments in this exercise are updated ver- 
sions of what is in the literature. Some arguments 
are new, as is the adoption of a single universal 
unit — the density parameter — that makes com- 
parisons across a broad variety of forms of energy 
immediate, but the central new development is 
that the considerable range of consistency checks 
demonstrates that many of the entries in the in- 
ventory are meaningful and believable. 

People have been making inventories for a long 
time. The medieval Domesday Book (1086-7) gave 
King William a picture of the wealth and organi- 
zation of the kingdom (and it gives us fascinat- 
ing insight into a society). The present-day cos- 
mic energy inventory similarly gives us a picture 
of the amount and organization of the material 
contents of the universe. It also offers us a way 
to assess the reliability of our picture, through 
checks of consistency. A early example of the lat- 
ter point is de Sitter's (1917) discussion of Ein- 



stein's (1917) static homogeneous universe. Under 
the assumption that observations can reach a fair 
sample of the universe, one can seek constraints 
on the cosmic mean mass density and space cur- 
vature and test the predicted relation between the 
two. Hubble's (1929) redshift-distance relation led 
to a revision of the predicted relation between the 
mass density and space curvature, and his use 
of redshifts to convert galaxy counts into num- 
ber densities greatly improved the estimate of the 
mean mass density (Hubble 1934).^ One may also 
consider the relations among the mean luminos- 
ity density of the galaxies, the production of the 
heavy elements, and the surface brightness of the 
night sky (Partridge & Peebles 1967; Peebles & 
Partridge 1967); the relation between galaxy col- 
ors, the initial mass function, and the star forma- 
tion history (Searle, Sargent, & Bagnuolo 1973); 
the relation between the last two sets of consid- 
erations (Tinsley 1973); the relation between the 
light element abundances and the baryon mass 
density (Gott et al. 1974); and the relation be- 
tween the luminosity density of the quasars and 
the mean mass density in quasars and their rem- 
nants (Sohan 1982). Basu and Lynden-Bell (1990) 
show how one can analyze what is learned from 
this rich set of considerations in terms of an en- 
tropy inventory. We have chosen instead to base 
this discussion on an energy inventory. 

Our inventory includes the mass densities in the 
various states of baryons. This is an updated ver- 
sion of the baryon budget of Fukugita, Hogan & 
Peebles (1998; FHP). Most entries in this part of 
the inventory have not changed much in the past 
half decade, while substantial advances in the ob- 
servational constraints have considerably reduced 
the uncertainties. It appears that most of the 
baryonic components arc observationally well con- 
strained, apart from the largest entry, for warm 
plasma, which still is driven by the need to bal- 
ance the budget rather than more directly by the 
observations. 

The largest entries, for dark matter and the cos- 
mological constant, or dark energy, are well con- 
strained within a cosmological theory that is rea- 



^Hubble's (1936) estimate based on galaxy masses derived 
from the velocity dispersion in the Virgo Cluster, which 
takes account of what is now termed nonbaryonic dark mat- 
ter, translates to density parameter Qm ~ 0.1, impressively 
close to the modern value. 



sonably well tested, but the physical natures of 
these entries remain quite hypothetical. We un- 
derstand the physical natures of magnetic fields 
and cosmic rays, but the theories of the evolution 
of these components, and the estimates of their 
contributions to the present energy inventory, are 
quite uncertain. The situation for most of the 
other entries tends to be between these extremes: 
the physical natures of the entries are adequately 
characterized, for the most part, and our estimates 
of their energy densities, while generally not very 
precise, seem to be meaningfully constrained by 
the observations. 

Several cautionary remarks are in order. First, 
some types of energy are not readily expressed as 
sums of simple components; we must adopt con- 
ventions. Second, there is no arrangement of cate- 
gories that offers a uniquely natural place for each 
component; we must again adopt conventions. 
Perhaps further advances in the understanding of 
cosmic evolution will lead to a more logically or- 
dered inventory. Third, it is arguably artificial 
to represent binding energies as very small nega- 
tive density parameters. The advantage is that it 
simplifies comparisons across the entire inventory. 
Fourth, the amount of space we devote is in accord 
with the importance of the issues of physics and 
astronomy, not with the value of the energy. For 
instance, the major components in the inventory 
are dark energy and dark matter, but they are 
physically simple in that their work is only grav- 
itational, so our discussion is rather short. Dissi- 
pativc energies are small in absolute values, but 
their physical significance is large, and our dis- 
cussions are on occasion rather lengthy. Finally, 
it is a task for future work to make some of our 
estimates more accurate by using data and com- 
putations that exist but are difficult to assemble. 
We mention the main examples in §3. 

The inventory, which is presented in Table 1, 
is arranged by categories and components within 
categories. The explanations of conventions and 
sources for each entry are presented in §2, in sub- 
sections with numbers keyed to the category num- 
bers in the first column of the table. Our dis- 
cussion of checks of the entries is not so simply 
ordered, because the checks depend on relations 
among considerations of entries scattered through 
the table. A guide to the considerable variety of 
checks detailed in §2 is presented in §3. 



Table 1 
The Cosmic Energy Inventory 



Components^ 



Totals" 



1 dark sector 

1.1 dark energy 

1.2 dark matter 

1.3 primeval gravitational waves 



0.954 ± 0.003 



0.72 ±0.03 
0.23 ±0.03 
< 10- 



-,-10 



2 primeval thermal remnants 

2.1 electromagnetic radiation 

2.2 neutrinos 

2.3 prestellar nuclear binding energy 



0.0010 ±0.0005 



10" 



3 baryon rest mass 

3.1 warm intergalactic plasma 

3.1a virialized regions of galaxies 

3.1b intergalactic 

3.2 intracluster plasma 

3.3 main sequence stars 
3.4 

3.5 white dwarfs 

3.6 neutron stars 

3.7 black holes 

3.8 substellar objects 

3.9 HI + Hel 

3.10 molecular gas 

3.11 planets 

3.12 condensed matter 

3.13 sequestered in massive black holes 



0.045 ± 0.003 



0.024 ± 0.005 
0.016 ±0.005 

spheroids and bulges 
disks and irregulars 



0.040 ± 0.003 



0.0018 ±0.0007 
0.0015 ± 0.0004 
0.00055 ±0.00014 
0.00036 ± 0.00008 
0.00005 ± 0.00002 
0.00007 ± 0.00002 
0.00014 ± 0.00007 
0.00062 ± 0.00010 
0.00016 ±0.00006 
10" 



-,-6 



10 
10" 



-5.6±0.3 



'(l + .„) 



4 primeval gravitational binding energy 

4.1 virialized halos of galaxies 

4.2 clusters 

4.3 large-scale structure 



-10" 
-10" 



5 binding energy from dissipative gravitational settling 

5.1 baryon-dominated parts of galaxies 

5.2 main sequence stars and substellar objects 

5.3 white dwarfs 

5.4 neutron stars 

5.5 stellar mass black holes 

5.6 galactic nuclei early type 

5.7 late type 



-10 

-10 
-10 
-10 
-10 
-10 
-10 



-8.8±0.3 
-8.1 

7.4 

5.2 

-5.8^ 



-10" 



6. poststellar nuclear binding energy 

6.1 main sequence stars and substellar objects 

6.2 diffuse material in galaxies 

6.3 white dwarfs 

6.4 clusters 

6.5 intergalactic 



_-|^Q-6.2±0.5 



-10" 



7 

7.1 

7.2 

7.3 

7.4 

7.5 

7.6 



poststellar radiation 

resolved radio-microwave 

far infrared 

optical 

X-7 ray 

gravitational radiation 



stellar mass binaries 
massive black holes 



-|^Q-10.3±0.3 

10-"! 

]^g-5.S±0.2 
]^Q-7.9±0.2 

^g-9±l 

1 |-j — 7.5ib0.5 



10" 



8 stellar neutrinos 

8.1 nuclear burning 

8.2 white dwarf formation 

8.3 core collapse 



-6.8 
-7.7 
-5.5 



10" 



cosmic rays and magnetic fields 



O + 0.6 
■^-0.3 



10 



kinetic energy in the intergalactic medium 



10" 



'^Bascd on Hubble parameter h — 0.7. 



2. The Energy Inventory 

The inventory in Table 1 assumes the now stan- 
dard relativistic Friedmann-Lemaitre ACDM cos- 
mology, in which space sections at fixed world 
time have negligibly small mean curvature,^ Ein- 
stein's cosmological term, A, is independent of 
time and position, the dark matter is an initially 
cold noninteracting gas, and primeval departures 
from homogeneity are adiabatic, Gaussian, and 
scale-invariant. Physics in the dark sector is not 
well constrained: A might be replaced with a dy- 
namical component,'^ as in the models for dark 
energy now under discussion, the physics of the 
dark matter may prove to be more complicated 
than that of a free coUisionless gas, and the initial 
conditions may not be adequately approximated 
by the present standard cosmology. If such com- 
plications were present we expect their effects on 
entries that are sensitive to the cosmological model 
would be slight, however, because the cosmological 
tests now offer close to compelling evidence that 
the ACDM model is a useful approximation to re- 
ality (Bennett et al. 2003a; Spergel et al. 2003; 
Tegmark et al. 2004a; and references therein). 

To help simplify the discussion we adopt a nom- 
inal distance scale, corresponding to Bubble's con- 
stant^ 



Ho == 70 km s"^ Mpc 



(1) 



The energy density, pi, in the form of component 
i is expressed as a density parameter. 



n. 



SnGp^ 
~3Hf' 



(2) 



Since Ho is at best measured to ten percent accu- 
racy (Freedman et al. 2001), an improved distance 
scale could produce noticeable revisions to the in- 
ventory. 



^The most compelling evidence is from the position of the 
first peak in the anisotropy power spectrum of the of cos- 
mic microwave background temperature field. The WMAP 
measurements (Bennett et al. 2003a) indicate Qm + ^x = 
1.02 ± 0.02. In this paper we simplify the discussion by 
assuming strictly flat space sections. 

•^The current limit on the index of the equation of state 
for the dark energy is ui = p/p < —0.78 at 95% (Bennett 
2003a); the bound w = — 1.02_Qj^g is obtained from the 
Type la supernova Hubble diagram under the assumption 
of flat space curvature (Riess et al. 2004) 

■^We also write Ho = lOOh km s~ Mpc~^, where conve- 
nient, but all entries in the inventory in Table 1 assume 
h = 0.7. 



The second column from the right in Table 1 
lists the density parameters in the components, 
and the last column presents the total for each 
category. Both columns sum to unity. The state- 
ment of errors requires some explanation. Most of 
the errors in categories 1 through 3 are rather well 
documented. The uncertainties in entries 2.1 and 
2.3 are too small to be relevant for the purpose of 
this inventory. Elsewhere a numerical value stated 
to one significant figure after the decimal might 
be supposed to be reliable to roughly ±0.1 dex, or 
about 30%. Where no digit is presented after the 
decimal point we hope our estimate might within 
a factor of ten of the real value. 

2.1. The Dark Sector 

The components in category 1 interact with the 
contents of the visible sector only by gravity, as far 
as is now known. This makes it difficult to check 
whether the dark energy — or Einstein's cosmo- 
logical constant, A — and the dark matter really 
have the simple properties assumed in the ACDM 
cosmology. Future versions of the inventory might 
contain separate entries for the potential, kinetic 
and gradient contributions to the dark energy den- 
sity, or a potential energy component in the dark 
matter. 

There is abundant evidence that the total mass 
density — excluding dark energy — is well be- 
low the Einstein-de Sitter value. That means, 
among other things, that the consistency of cross- 
checks from the many ways to estimate the mass 
density provides close to compelling evidence that 
the gravitational interaction of matter at distances 
up to the large-scale flows is well approximated 
by the inverse square law, and that starlight is 
a good tracer of the mass distribution on scales 
> 100 kpc.5 

An example that illustrates the situation, and 
will be used later, is the estimate from weak lens- 
ing of the mean galaxy surface mass density con- 
trast, 

S(<y)-E(y)=A(%/lMpc)-", (3) 



®For early discussions see Peebles 1986 and Bahcall, Lubin 
& Dorman (1995). Recent reviews of the situation are in 
Pukugita (2001), Peebles & Ratra (2003), and Bennett et 
al. (2003a). The present observational situation is reviewed 
in footnote 6. 



where S(y) is the ensemble average surface mass 
density at projected distance y from a galaxy and 
E(< y) is the mean surface density within distance 
y. The measurements by McKay et al. (2001) 
yield A = 2.5tolhmQ pc'^ and a = -0.8 ± 0.2, 
and they indicate that the power law is a good 
approximation to the measurements in the range 
of projected radii 

70 kpc < y < 1 Mpc. (4) 

If the galaxy autocorrelation function, 

^W = (^)'' , 7 = 1.77, ro = 5/1-1 ^jp^^ (5) 

is a good approximation to the galaxy-mass cross 
correlation function on the range of scales in equa- 
tion (4), then the mean surface density is 



^(y) 



H-y 



r(i/2)r((7-i)/2) 
r(7/2) 



(6) 



This agrees with equation (3) if the density pa- 
rameter belonging to the mean density pm of the 
mass that clusters with the galaxies is 



J7m(weak lensing) — 0.20^ 



0.05- 



(7) 



This is in the range of estimates of the value of this 
parameter now under discussion, consistent with 
the assumption that galaxies are useful tracers of 
mass. 

Our adopted value for the total mass density in 
nonrelativistic matter, dark plus baryonic, is 



u ^YY) " ^ 



DM 



Vt^ + n^ = 0.28 ±0.03. 



(8) 



This is in the range of most current estimates.^ 
The measurement may not be tightly constrained. 



"Spergel et al. derived Q.mh? = 0.13 - 0.14 from WMAP 
either with or without using the constraint from the power 
spectrum of the 2dF galaxy distribution. This gives f2m = 
0.265-0.286 at our fiducial h = 0.7. Tegmark et al. (2004a) 
obtained the central value ^rnh? = 0.14 with WMAP 
data alone, and 0.145 with the constraint from the SDSS 
galaxy clustering. For h = 0.7 these values are respec- 
tively Um = 0.29 and 0.30. The estimate from the 2dF- 
GRS power spectrum (W. Percival and the 2DFGRS team, 
2004, private communication) is Qm,h = 0.164 ± 0.016, or 
^m = 0.23 ± 0.02 at our distance scale. This is 1.4 stan- 
dard deviations below equation (8). The lower estimate of 
Clm is due to the somewhat larger value of the small-scale 
power spectrum compared to SDSS. Note that departures 



however, and there is the possibility of adjustment 
of this important parameter beyond the error flag 
in equation (8). 

We make use of the fact that equation (3) is 
close to the limiting isothermal sphere mass dis- 
tribution, 

p{r) = a^/2nGr'^ (9) 

If we connect this form to the power law in equa- 
tion (5) at the nominal virial radius ry defined by 
p{< ry)/pc = 200, we obtain 



r„ = 220h^\pc, 



a — 160 km s 



(10) 



This measure of the characteristic one-dimensional 
velocity dispersion, cr, in luminous galaxies agrees 
with the mean in quadrature. 



(cr2)i/2 ^ 160 km s-\ 



(11) 



weighted by the FHP morphology fractions, of the 
dispersions a = 225 km s""'^ for elliptical galaxies, 
a = 206 km s^^ for SO galaxies (de Vaucouleurs & 
Olson 1982) and a — 136 km s"^ for spiral galax- 
ies (Sakai et al. 2000), all at luminosity Lb = L*^. 
The isothermal sphere model defines a character- 
istic mass My within r^. The ratio of My to the 
characteristic galaxy luminosity (eq. [(53] below) 



IS 



MyjLr = 180/1, My I Lb = 250/i, 



(12) 



in Solar units. This is consistent with the esti- 
mate of MjL within 220h^^ kpc from weak lens- 
ing shear (McKay et al. 2001). '^ 

The product of My with the effective number 
density of luminous (L*) galaxies, Ug — 0.017/1"^ 
Mpc""^ (eq [52]), gives an estimate that the mean 



from the standard assumptions, including fiat space curva- 
ture, scale-invariant initial conditions, and neglible tensor 
perturbations, could lead to changes beyond the quoted er- 
rors. We also refer to estimates from the Type la sypernova 
redshift-magnitude relation, flm = 0.29_gg3 (Riess et al. 
2004), and from dynamics, including Clm = 0.17it0.05 from 
the cluster mass function as a function of redshift (Bahcall 
& Bode 2003) and f^^ = 0.30 ±0.08 from the redshift space 
two-point correlation function (Hawkins et al. 2003), under 
the assumption that the bias parameter is fe = 1 (Verde et 
al. 2002). 
^It is worth noting that equation (12) is not far from 
Zwicky's (1933) dynamical estimate for the Coma Cluster, 
M/L ^ 100 at our distance scale. 



mass fraction within the virial radii of the large 
galaxies, 

^il^ = 0.6. (13) 

That is, we estimate that about 60% of the dark 
matter is gathered within the virialized parts of 
normal galaxies. 

We consider now how the estimate of ilm in 
equation (8) compares to the estimate from the 
mass-to-light ratio and the integrated mean lu- 
minosity density. The estimates from the SDSS 
broad-band galaxy luminosity functions are (Ya- 
suda et al. in preparation; see also Blanton et al. 
2003) 

Cb = (1.9±0.2) X lO*/^L0Mpc"^ 
Cr = (2.3±0.2) X lO^/^L0Mpc^ (14) 
A = (3.6±0.4) X lO^/iLoMpc^^ 

These densities are the values at z w 0.05. The 
value of Cb is inferred from the densities in other 
color bands. ^ The luminosity density in the z band 
is quoted for the later use. The product of the lu- 
minosity density with M / L in equation (12) yields 
the density parameter in matter within the virial 
radii of galaxies, ^m.v ~ 0.18 and 0.16 for the B 
and r bands, respectively. On dividing by equation 
(13) we arrive at Vl^ = 0.31 and 0.27, consistent 
with equation (8). 

Note that M/ L in equation (12) is about half 
the value estimated for rich clusters, (M/L)ciuster = 
450 ± 100 (FHP). In rich clusters all mass in the 
outskirts of galaxies, beyond their virial radii, is 
integrated in the mass estimate, so that it is rea- 
sonable to suppose that it gives a larger value by 
the inverse factor in equation (13). When the 
cluster value for M/L is multiplied by the B band 
luminosity density we get Q.m = 0.32, consistent 
with the other estimates. 

Entry 1.3 assumes inflation has produced grav- 
itational waves with a scale-invariant spectrum, 
meaning the strain 5 appearing at the Hubble 
length is independent of epoch. The density pa- 
rameter associated with gravitational waves with 
wavelengths on the order of the Hubble length is 



Vlg ^ (5^, and the absence of an appreciable ef- 
fect of gravitational waves on the anisotropy of 
the 3 K thermal cosmic background radiation in- 
dicates 5 < 10^'"'. The gravitational waves pro- 
duced by cosmic phase transitions, if detected or 
convincingly predicted, might be entered in this 
category. Gravitational waves from the relativis- 
tic collapse of stars and galactic nuclei are included 
in category 7. 

The other entries in the first category in Table 1 
are computed from equation (8) and our estimates 
of the other significant contributions to the total 
mass density, under the assumption that the den- 
sity parameters sum to unity, that is, space curva- 
ture is neglected. 

2.2. Thermal Remnants 

2.2.1. Cosmic Background Radiation 

Entry 2.1 is based on the COBE measurement 
of the temperature of the thermal cosmic elec- 
tromagnetic background radiation (the CMBR), 
To = 2.725 K (Mather et al. 1999). The COBE 
and UBC measurements (Mather et al. 1990; 
Gush, Halpern, & Wishnow 1990) show that the 
spectrum is very close to thermal. It has been 
slightly disturbed by the observed interaction with 
the hot plasma in clusters of galaxies (LaRoque et 
al. 2003 and references therein). The limit on 
the resulting fractional increase in the radiation 
energy density is (Fixsen et al. 1996) 



6u/u = 4j/ < 6 X 10 



-5 



(15) 



'The ratios of luminosity densities, 1 : 1.20 : 1.87, are 
approximately what is expected from the average colors, 
{B - r) = 1.0 and {r - z) = 0.6 [{B - r)Q = 0.82, 
(r-^)0 = O.12]. 



This means that the background radiation en- 
ergy density has been perturbed by the amount 
Afl < 10"^'^. Improvements of this number are 
under discussion (e.g. Zhang, Pen & Trac 2004), 
and might be entered in a future version of the 
inventory. 

The thermal background radiation has been 
perturbed also by the dissipation of the primeval 
fluctuations in the distributions of baryons and ra- 
diation on scales smaller than the Hubble length 
at the epoch of decoupling of baryonic matter and 
radiation. If the initial mass fluctuations are adi- 
abatic and scale-invariant the fractional perturba- 
tion to the radiation energy per logarithmic incre- 
ment of the comoving length scale is Su/u ~ (5^, 
where Sh ^ 10^^ is the density contrast appear- 
ing at the Hubble length. This is small compared 



to the subsequent perturbation by hot plasma 
(eq [15]). 

Entry 2.2 uses the standard estimates of 
the rehct thermal neutrino temperature, T^ — 
{4:/ll)^'^To, and the number density per family, 
rii, = 112 cm~^. We adopt the neutrino mass dif- 
ferences from oscillation experiments (Fukuda et 
al. 1998; Kameda et al. 2001; Eguehi, et al. 2003; 
BahcaU & Pefia-Garay 2003), 



ml^ ~ ml^ = 0.002 ± 0.001 eV^ 
ml^ - ml^ = 6.9 x 10"^ eV^ 



(16) 



where the neutrino mass eigenstates are ordered as 
mjyj < TOjyj < 771^3. Entry 2.2, the density param- 
eter rii, in primeval neutrinos, assumes that m^^ 
may be neglected. The upper limit from WMAP 
and SDSS is fl^ < 0.04 (Tegmark et al. 2004a). 
At this limit the three families would have almost 
equal masses, m^ = 0.6 eV. This may not be very 
likely, but one certainly must bear in mind the 
possibility that our entry is a considerable under- 
estimate. 

2.2.2. Primordial Nucleosynthesis 

Light elements are produced as the universe ex- 
panded and cooled through kT ~ 0.1 MeV, in 
amounts that depend on the baryon abundance. 
The general agreement of the baryon abundance 
inferred in this way with that derived from the 
CMBR temperature anisotropy gives confidence 
that the total amount of baryons — excluding 
what might have been trapped in the dark mat- 
ter prior to light element nucleosynthesis — is se- 
curely constrained. 

Estimates of the baryon density parameter from 
the WMAP and SDSS data (Spergel et al. 2003; 
Tegmark et al. 2004a), and from the deuterium 
(Kirkman et al. 2003) and helium abundance 
measurements (Izotov & Thuan 2004) are, respec- 
tively, Ob/i^ == 0.023 ± 0.001, 0.0214 ± 0.0020, 
and 0.0131;q°q§]^, where the last number is the all- 
sample average for helium from Izotov & Thuan. 
We adopt 



n^h^ =0.0225 ±0.0015, 



(17) 



close to the mean of the first two. Since the 
relation between the helium abundance and the 
baryon density parameter has a very shallow slope. 



an accurate abundance estimate (say, with < 1% 
error) is needed for a strong constraint on fl\yh^. 
We consider that the current estimates may still 
suffer from systematic errors which are not in- 
cluded in the error estimates in the literature.^ 
Within the standard cosmology our adopted value 
in equation (17) requires that the primeval helium 
abundance is 



r„ = 0.248 ±0.001, 



(18) 



and the ratio of the total matter density to the 
baryon component is 



r2„/f7t = 6.11 ±0.89. 



(19) 



We need in later sections the stellar helium pro- 
duction rate with respect to that of the heavy ele- 
ments. The all-sample analysis of Izotov & Thuan 
(2004) gives AY/AZ ~ 2.8 ± 0.5. The value 
derived by Peimbert, Peimbert & Ruiz (2000) 
corresponds to 2.3 ± 0.6. These value may be 
compared to estimates from the perturbative ef- 
fects on the effective temperature-luminosity rela- 
tion for the atmosphere of main sequence dwarfs, 
AY/AZ ~ 3 ± 2 (Pagel & Portinari 1998), and 
2.1 ± 0.4 (Jimenez et al. 2003). From the initial 
elemental abundance estimate in the standard so- 
lar model of BahcaU, Pinsonneault & Basu (2001; 
hereinafter BP2000) we derive AY/AZ ~ 1.4. We 
adopt 

Ar/AZ~2±l. (20) 

Nuclear binding energy was released during nu- 
cleosynthesis. This appears in entry 2.3 as a neg- 
ative value, meaning the comoving baryon mass 
density has been reduced and the energy density 
in radiation and neutrinos increased. The effect on 
the radiation background has long since been ther- 
malized, of course, but the entry is worth recording 



^We note, as an indication of the difficulty of those obser- 
vations, that helium abundances inferred from the triplet 
4d-2p transition (A4471) are lower than what is indicated 
by the triplet 3d-2p (A5876) and singlet 3d-2p (A6678) 
transitions, by an amount that is significantly larger than 
the quoted errors. Another uncertainty arises from stellar 
absorption corrections, which are calculated only for the 
A4471 line. The table given in Izotov and Thuan suggests 
that a small change in absorption corrections for the A4471 
line induces a sizable change in the final helium abun- 
dance estimate. We must remember also that the value 
of AY/AZ, which is needed to derive Yp, is not very well 
determined. 



for comparison to the nuclear binding energy re- 
leased in stellar evolution. For the same reason, we 
compute the binding energy relative to free pro- 
tons and electrons. The convention is artificial, 
because light element formation at high rcdshifts 
was dominated by radiative exchanges of neutrons, 
protons and atomic nuclei, and the abundance of 
the neutrons was determined by energy exchanges 
with the cosmic neutrino background. It facili- 
tates comparison with category 6, however. The 
nuclear binding energy in entry 2.3 is the product 



Mz — M* -I- 3 for a galaxy sample at z « 0.05. 
Our estimate of the resulting luminosity-function- 
weighted mean is 



-riNB,Hc = 0.0071 ypf}b = 10 



-4.1 



(21) 



This is larger in magnitude than the energy in the 
CMBR today. 

2.3. The Baryon Rest Mass Budget 

The entries in this category refer to the baryon 
rest mass: one must add the negative binding en- 
ergies to get the present mass density in baryons. 
The binding energies are small and the distinction 
purely formal at the accuracy we can hope for in 
cosmology, of course, with the conceivable excep- 
tion of the baryons sequestered in black holes. 

We begin with the best-characterized com- 
ponents, the stars, star remnants, and planets. 
We then consider the diffuse components, and 
conclude this subsection with discussions of the 
baryons in groups and the intergalactic medium 
and the lost baryons in black holes. 

2.3.1. Stars 

This is an update of the analysis in FHP. Fol- 
lowing the same methods, we estimate the baryon 
mass in stars from the galaxy luminosity density 
and the stellar mass-to-light ratio, Mstars/L, along 
with a stellar initial mass function (the IMF) that 
allows us to estimate the mass fractions in various 
forms of stars and star remnants. 

Kauffmann et al. (2003) present an extensive 
analysis of the stellar mass-to-light ratio based on 
ugriz photometry for 10^ SDSS galaxies and a 
population synthesis model (Bruzual & Chariot 
2003) that is meant to take account of the stel- 
lar mctallicitics and the star formation histories. 
Their estimate of the stellar mass-to-light ratio is 
Msturs/Lz — 1-85 for luminous galaxies, with mag- 
nitudes Mz < M* — 0.8, and it decreases gradu- 
ally to Mstars/Lz — 0.65 for fainter galaxies with 



(Mstars/i.) ^ 1.5 ± 0.3, 



(22) 



for the IMF Kauffmann et al. used. Equation (22) 
represents the present mass in stars and stellar 
remnants, and it excludes the mass shed by evolv- 
ing stars and returned to diffuse components. 

The estimate of Msturs/L assumes a universal 
IMF, which is not thought to seriously violate the 
observational constraints. We note, however, that 
a possible change of the IMF at very high redshift 
need not seriously affect our analysis because star 
formation a high redshift contributes little to the 
present total mass in stars. The IMF is particu- 
larly uncertain at the subsolar masses that make 
little contribution to the light but can make a con- 
siderable contribution to the mass. We consider 
two continuous broken power law models, of the 
form 

diV/dm oc m"("'+i\ (23) 

where, in the first model, 

X = -0.5, 0.01 <m < O.lniQ, 

= 0.25 0.1 < m < IniQ, (24) 

= 1.35 1< m < lOOm0, 

and, in the second model, 

X = -0.7 0.01 <m< O.OSmg, 

= 0.3 0.08 <m< O.5m0, (25) 

= 1.3 0.5 < m < lOOm©. 

The first line in the first model is from Burgasser 
et al. (2003), the second line is from Rcid, Gizis, & 
Hawley (2002), and the third line, for m > Im©, 
is the standard Salpeter (1955) IMF. The second 
model is from Kroupa (2001). Yet another IMF, 
that of Chabrier (2003), is in between these two. 
For our model IMF we take the Salpeter index 
for TO > Itoq. It is known that the Salpeter 
slope gives satisfactory UBV colors and Ha equiv- 
alent widths for normal galaxies (Kennicutt 1983), 
whereas an IMF with a steeper slope (e.g., Scalo 
1986) is not favored in this regard. The consensus 
seems to be that the power law index is smaller 
than the Salpeter value at sub-Solar masses. The 
two IMF given above still differ significantly at 



m < niQ, however, For the subsolar IMF, we take 
the geometric mean of the above two models after 
mass integration, and we take the difference as an 
indication of the error (±18%). 

The stehar mass-to-light ratio in equation (22) 
assumes the IMF in equation (25). With our 
adopted IMF the stellar mass-to-light ratio is 1.18 
times the number in equation (22).-'^'' Thus we get 
our fiducial estimate, 



(Mstar/i.) = 1.23 ±0.33. 



(26) 



This translates to M/Lb ~ 2.4, or 0.7 times that 
used in FHP, which employed the subsolar mass 
IMF of Gould, Bahcall & Flynn (1996). In the 
Salpeter IMF, with x = 1.35 cut off at O.lm©, 
the mass-to-light ratio is 1.48 times our adopted 
value. The Kennicutt (1983) IMF resuhs in 0.81 
times equation (26). 

The IMF at substellar masses, m < O.O8m0, 
is more uncertain, but recent observations of T 
dwarfs in the solar neighbourhood indicate x < 
(e.g., Burgasser et al. 2003; 2004). In the two 
IMF models quoted above the substellar mass is 6 
to 9% of the mass integral for 0.08 < m < Ittiq. 
We adopt 8% and assign an error of 50%. 

We estimate the mass density locked in stars, 
including those in dead stars, to be Ogtars — 
0.0024 ±0.0007. A comparable estimate is derived 
from bj, J, Kg multicolour photometry of 2MASS 
combined with 2dF data by Cole et al. (2001), 
figta^^g = 0.0029 ± 0.0004 with the IMF we have 
adopted. For the energy inventory we take the 
mean of our present number and that of Cole et 
al.: 

r^.tars = 0.0027 ± 0.0005. (27) 

This means that the stars contain 6.0 ± 1.3% of 
the total baryons. The FHP estimate is figtars — 
0.0019 - 0.0057. Equation (27) also is consis- 
tent within the errors with the more recent es- 
timates by Salucci & Persic (1999), Kochanek et 
al. (2001), and Glazebrook et al. (2003), and with 
ShuU's (2003) baryon inventory. 

We attempt to partition the stars into their 
species. Our estimates of the mass fractions in 
stars on the main sequence (MS) and substellar 



^"The IMF are normalised so that the mass integrals between 
0.9 and 2.0mQ arc equal. The result is virtually identical 
with those with the two IMFs normalised at IrrtQ . 



objects (SS), and the mass fractions in stellar rem- 
nants, including white dwarfs (WD), neutron stars 
(NS) and stellar mass black holes (BH), are shown 
in Table 2. Stars on the MS are represented by 
the present-day mass function (PDMF) , which for 
1 < TO < 100to,q is constructed by multiplying the 
IMF by m^^-^ to account for the lifetime of mas- 
sive main sequence stars; for m < Ito-0 the PDMF 
is the same as the IMF. A more detailed estimate 
of the relation between PDMF and IMF is possi- 
ble but not needed for our purposes because the 
mass fraction in stars with masses to. > TOq is 
small. We take to = O.O8to0 as the mass divid- 
ing hydrogen-burning and substellar mass 'brown 
dwarfs' (Hayashi & Nakano 1963; Kumar 1963) 
and the observational limit O.O1to,0 as the lower 
end of substellar masses. 

As indicated in the third line of Table 2, we 
take the average mass of a white dwarf to be 
(to,) = 0.62 ?Ti0, for consistency with the model we 
describe in section 2.5.2 (eq. [70]). This is close 
to the value from recent observations, 0.604 ?Ti0 
(Bergeron & Holberg, in preparation; Bergeron 
private communication) . We assume all stars with 
initial masses in the range 1 < tti < 8 to,0 end up 
as white dwarfs with the adopted average mass, 
and the rest of the mass returns to the interstel- 
lar medium. With our adopted IMF this model 
predicts that the ratio of masses in white dwarfs 
to main sequence stars is piy^TD) / piMS) — 0.176. 
This can be compared to the observations. The 
2dF survey (Vennes et al. 2002), with the 
use of our mean mass 0.62 ?Ti0, yields a mea- 
sure of the mass density of DA white dwarfs, 
/9(DA-WD,local) = (4.2 ± 2.3) x \Q-^mQ])C-'^ . 
This is multiplied by 1.3 to account for DB, DQ 
and DZ white dwarfs (Harris et al. 2003), to give 
the total white dwarf mass density p(WD, local) = 
(5.5±3.0) X 1O-3to0 pc-3. This divided by the lo- 
cal density of main sequence stars, /9(MS, local) = 
0.031 ± 0.002 TO0 pc-3^ from Reid, Gizis, & 
Hawley (2002), and 0.041 ± 0.003 TO0 pc'^ from 
Chabrier (2003), gives p{YfV)) / p{MS) = 0.16 ± 
0.10, which agrees with our model prediction, al- 
beit with a large uncertainty. 

We assume all stars in the initial mass range 



8 < TO < 25 TO, 







(28) 



end up as neutron stars with mass 1.35 to,0 (Nice, 
Splaver, & Stairs 2003; Thorsctt & Chakrabarty 



Table 2 
Star mass fractions 



initial 



mass range 



fate 



fraction 



mass consumed 



0.01 < m < 0.08 


SS 




0.052 


0.052 


0.08 < m < lOO"" 


MS 




0.769 


0.769 


1 < m < 8 


WD 


0.62 


0.135 


0.463 


8 < m < 25 


NS 


1.35 


0.019 


0.186 


25 < m < 100 


BH 


7.5 


0.025 


0.146 



1.0 



1.616 



^Mass in units of Solar masses 

The gas consumed to make unit mass of stars now present. 
^For the mass range 1 < m < IOOtmq those stars burning today are counted. 



1999), and the rest of the star mass is recycled. Es- 
timates of the lower critical mass for stellar core 
collapse vary from 6 to 10 mQ (Reimers & Koester 
1982; Nomoto 1984; Chiosi, Bertelh, & Bressan 
1992). The upper critical mass for the formation 
of a neutron star is more uncertain; our choice in 
equation (28) is taken from Heger et al. (2003). 
We assume all stars with main sequence mass 
above the limit in equation (28) end up as stellar 
black hole remnants with mass m, ~ 7.5 toq, with 
the rest of the mass recycled. This again follows 
Heger et al. The remnant mass is quite uncertain, 
and Heger et al. also indicate that some stars with 
masses above the range in equation (28) may pro- 
duce neutron stars. 

The last column in Table 2 is the cumulative 
amount of mass that has gone into stars, including 
what is later shed, normalised to the mass in col- 
umn (4). According to these estimates, 1.6 times 
the mass in observed stars was used (and reused) 
in star formation. 

Entries 3.3 to 3.8 in Table 1 are based on the 
partition of fistars among main sequence stars and 
star remnants in the fourth column of Table 2, 
with the FHP partition between the two galaxy 
population types in the proportion 



spheroid : disc — 0.74 : 0.26 



(29) 



of mass in ellipticals plus bulges of spirals to the 
mass in disks plus irregulars. 

2.3.2. Consistency with the Star Formation Rate 

We can compare our estimate of the mass in 
stars to the mass integrated over the star forma- 



tion rate density (the SFR) from high redshift to 
the present. The Ha luminosity density at zero 
redshift measured from an SDSS nearby galaxy 
sample is C{YLa) = 10^^-^^-o:Js /j erg s^^Mpc^^ 
(Nakamura et al. 2004). This agrees with the 
earlier value obtained by Gallego et al. (1995). 
Within the star formation models Glazebrook et 
al. (1999) explored, Im© yr^^ of star formation 



+0.92^ 
-0.24; 



X lO'' 



produces Ha luminosity Lq — (2.00 

erg s^^ for our adopted IMF. The ratio of £(Ha) 

to La is an estimate of the present-day SFR, 



V'(O) - 0.007lti;:^^^^mo yr-iMpc-\ 



(30) 



The evolution of the SFR at z < 1 is now ob- 
servationally well determined (Lilly et al. 1996; 
Heavens et al. 2004; see Glazebrook et al. 1999 
for a summary of the measurements up to that 
time), and may be approximated as 



V'(t) = V(0)e 



(to-t)/r,f 



z < 0.85, 



where 



Tsf ~ 2.9 Gyr, 



(31) 



(32) 



and to ^ ^ is the time measured back from the 
present in the ACDM cosmology (Fukugita & 
Kawasaki 2003; see also Glazebrook et al. 1999). 
The situation at higher redshift is still controver- 
sial (Madau et al. 1996; Steidel et al. 1999). 
We shall suppose that the SFR is constant from 
z = 0.85 back to the start of star formation at 
redshift Zf, at the value 

V- = 0.0877 TO0 yr"lMpc"^ 0.85 < z < zj. (33) 
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In this model the integrated star formation is 



^s 



0.0038 if Zf == 2, 
0.0047 if Zf = 3, 
0.0053 if Zf = 5. 



(34) 



There is a significant downward uncertainty, 
(Jfigtars ~ —0.0025, arising from the uncertainty 
in the star formation rate (eq.[30]). Equation (27) 
and Table 2 indicate that the mass processed 
into stars, including what was later shed, is 
n = 0.0027 X 1.62 = 0.0044. This is consis- 
tent with the cumulative star formation in equa- 
tion (34). There would be a problem with the 
numbers if there were reason to believe that the 
SFR continued to increase with increasing redshift 
well beyond z = 1. A constant or declining SFR 
at 2 > 1 is well accommodated with our estimates 
for the baryon budget and current ideas about the 
star formation history at high redshift. 

The star formation history determines the av- 
erage reciprocal redshift factor (as in eq. [102]) 
for the effect of redshift on the integrated comov- 
ing energy density of electromagnetic radiation 
(and other forms of relativistic mass) generated 
by stars. In our model for the SFR the correction 
factor is 



((1 + ^)-')- 



^dtC 



JdtC/{l + z) 
2.0 ±0.15, 



(35) 



where C{t) is the bolomctric luminosity density 
per comoving volume, which we are assuming is 
proportional to the star formation rate density. 
The numerical value assumes the redshift cutoff 
is in the range Zf ~ 2 to 5. For the purpose of 
comparison of the accumulation of stellar products 
to the present rate of production, another useful 
quantity represents the integrated comoving den- 
sity of energy radiated in terms of the effective 
time span normalised by the present-day luminos- 
ity density. 



Aieff == Cita)-^ I dtC{t) = 85 ± 13 Gyr, 



(36) 



The value assume the range of models for the star 
formation history in equation (34). 

We can compare the rate of stellar core collapse 
in our model to observations of the supernova rate. 



(See Fukugita & Kawasaki 2003 and Madau, della 
Valle & Panagia 1998 for similar analyses.) The 
present SFR in equation (30) and the critical min- 
imum mass 8 m© in equation (28) imply that the 
present rate of formation of neutron stars and stel- 
lar black holes is expected to be 



R = 



Jg dm dN /dm 
Jp gg dm m dN/dm, 
0m79troT9 (lOOyr)-iMpc 



(37) 



Our estimate of the observed supernova rate is 

i?SN = Omref'oZlt (100yr)-iMpc-3. (38) 

This is the geometrical mean of the rates from 
three surveys for Type II and Type Ib/c su- 
pernovae, 0.037, 0.018 and 0.017, in units of 
/i^(100yr)~iMpc"^ (Tammann, LofHer & Schroder 
1994; CappeUaro et al. 1997; van den Bergh & Mc- 
Clure 1994; see also CappeUaro, Evans & Turatto 
1999). We conclude that, if most stars with initial 
masses greater than about 8 ttiq produced Type 
II and Ib/c supernovae, then our model for the 
star formation history would pass this consistency 
check. 

We remark that in our model the comoving 
number density of supernovae of types II and Ib/c 
integrated back to the start of star formation is 



i?Ateff- = 7 ± 3 X W Mpc"^ 



(39) 



where the effective time span is given in equation 
(36). 

2.3.3. Planets and Condensed Matter 

The mass in planets that are gravitationally 
bound to stars must be small, but it is of partic- 
ular interest to us as residents of a planet. Marcy 
(private communication; see also Marcy & Butler 
2000) finds that about 6.5% of nearby FGKM stars 
have detected Jovian-like planets, and that an ex- 
trapolation to planets at larger orbital radii might 
be expected to roughly double this number. In 
our model for the PDMF the ratio of the number 
density of stars in the mass range 0.08 to 1.6 m© 
to the mass density in stars is n/ p ~ 2.1mQ^. The 
product of this quantity with the mass density in 
stars (eq. [27]), the fraction 0.13, and the ratio of 
the mass of Jupiter to the Solar mass is 



n 



planets 



10 



-6.1 



(40) 
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Marcy indicates that stars with lower metalhcity 
have fewer planets, but that may not introduce a 
serious error because there are fewer low metalhc- 
ity stars. 

There is a population of interstellar planets that 
have escaped from stars that have rapidly shed 
considerable mass. If the number of planets per 
star is independent of the stellar mass then, in 
our IMF, stars in the mass range 1 to 1.6toq had 
0.025 times the number of stars that are bound to 
main sequence stars. Among them about 2/3 are 
likely to be swallowed by the host stars in their gi- 
ant star phases. This leaves the mass in liberated 
planets, and those still bound to white dwarfs, 
at about 0.07 times the mass density in equa- 
tion (40). Planets associated with more massive 
stars may more likely escape during core collapse, 
but the numbers of stars, and perhaps planets, is 
small. The large uncertainty in equation (40) is 
whether the stellar neighborhood is a fair sample. 
This leads us to enter the order of magnitude in 
Table 1. 

We consider separately the mass in objects 
small enough to be held together by molecular 
binding energy rather than gravity. The dom- 
inant amount of the former is interstellar dust. 
It is known (Drainc 2003) that > 90% of silicon 
atoms in interstellar matter are condensed into 
grains, probably dominantly in the forms of cn- 
statite (MgSiOs) and forsterite (Mg2Si04) with 
the former likely three to four times more abun- 
dant (Molster et al. 2002). This means the mass in 
this form is Z (silicates) /Z — 0.17 times the mass 
in heavy elements. Drainc (2003) suspects that the 
dominant contribution to the carbonaceous mate- 
rial is in the form of polycyclic aromatic hydrocar- 
bon with the inferred abundance C/H=6xl0^'''. 
Adding this to the silicates, we find that the mass 
fraction becomes Z(dust)/Z — 0.20. The for- 
mation of dust with iron or other forms of car- 
bonates could increase this number. Entry 3.12, 
^dust = 10^^'®, is the product of the density pa- 
rameter in cool gas (entries 3.9 plus 3.10) with the 
mean metalhcity discussed below (and displayed 
as eq. [88]) and the 20 percent heavy element mass 
fraction in dust. Our estimate, which assumes 
that what we know about the Milky Way applies 
to other galaxies, is crude, but it seems likely that 
the mass in dust exceeds the mass in planets. 

There are larger objects in which molecular 



binding is important. The gravitational binding 
energy of a roughly spherical object with mass m 
and radius r is ~ Gm? /r, and the molecular bind- 
ing energy is roughly 1 eV per atom. The molecu- 
lar binding energy is larger than the gravitational 
energy when 

m< 1.2x10^® g{^gcmr^/pf/'^/A, (41) 

where the mass density is p and the mean atomic 
weight is A. For silicates, this bound is ?Ti = 
8 X 10^^ g. That is, gravitational binding energy 
dominates in the Earth and molecular binding en- 
ergy dominates in the Moon and asteroids. Tru- 
jillo, Jcwitt, & Luu (2001) estimate that the mass 
in the Kuiper Belt objects is about a tenth of an 
Earth mass, or IO^^'^toq. A comparable mass is 
in the moons in the Solar system. If this mass 
fraction were common to all stars the density pa- 
rameter in these objects would be 

^asteroid ^ 10~ , (42) 

a small fraction of the mass in dust. 

2.S4. Neutral Gas 

The recent blind HI surveys are a significant 
advance over the data used by FHP to estimate the 
mass density in neutral atomic gas. The largest 
survey, HIPASS, with 1000 galaxies (Zwaan et al. 
2003), gives 



n 



HI 



4.2 ± 0.7 X 10-"*. 



(43) 



The increase over the value quoted in FHP (1.5 
times the upper end value of FHP) illustrates the 
advantage of blind surveys over observations of 
programmed galaxies. The molecular hydrogen 
abundance from the CO survey of Kcrcs, Yun & 
Young (2003) is 



n^ 



1.6 ± 0.6 X 10" 



(44) 



The sum of these two values is multiplied by 1.38 
to accommodate helium. 

We place in entry 3.9 the atomic hydrogen and 
the helium abundance belonging to atomic and 
molecular hydrogen. Entry 3.10 is the molecular 
hydrogen component. The mass in this neutral 
gas is 1.7±0.4% of the total baryon mass. 
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2.3.5. Intracluster Plasma 

The estimate of the plasma mass in rich clus- 
ters of galaxies depends on a convention for the 
cluster radii and masses. We use the mass M200 
contained by the nominal virial radius, r2oo- This 
definition is more faithful to a physical definition 
of the part of a cluster that is close to dynamical 
equilibrium, and it also traces the X-ray radius 
which is the definition of our hot gas. In the lim- 
iting isothermal sphere model the relation to the 
mass within the AbcU radius ta = 1.5/i^^ Mpc is 

Ma = Ml/^Mll^, Mo = 1.1 X IQi^mQ. (45) 

This agrees with the estimates of the two masses 
given by Reiprich & Bohringer (2002), Table 4. 
We adopt the Reiprich & Bohringer cluster mass 
density parameter. 



O, 



0.012 



+0.03 
-0.04' 



(46) 



for the cluster mass limit 



M200 > 5 X IO^^toq, Ma > 1.3 x lO^^mQ. (47) 

That is, 4% of the mass is assembled in rich clus- 
ters. The Bahcall & Cen (1993) Abell mass func- 
tion is consistent with the Reiprich & Bohringer 
estimate. We caution that the Bahcall et al. 
(2003) SDSS mass function is about half what we 
adopt, perhaps because SDSS samples a subset of 
the clusters. Note also that the integral over the 
cluster mass within the Abell radius gives a sub- 
stantially larger value for fl^ because Ma > Mo 
at the masses M > lO^^r7i0 that dominate the 
integral for a given mass function. 

The abundance of hot baryons in clusters is ob- 
tained from the expression 

^b,cl — 0.012(fib — ^stars)/^m 

= 0.0018 ±0.0007. (48) 

The change from the value in FHP is mainly due 
to the definition of the cluster mass, as also noted 
by Reiprich & Bohringer. 

2.3.6. Massive Black holes 

We follow the standard idea that the massive 
objects in the centers of galaxies are black holes 
that formed by the accretion of baryons. Baryons 
entering black holes are said to lose their identity, 



but for our accounting it is appropriate to consider 
them to be sequestered baryons. 

We define the characteristic efficiency factor e„ 
for black hole formation (where the subscript is 
meant to distinguish the massive black holes in the 
centers of galaxies from stellar mass black holes) 
by the black hole mass produced out of an initial 
diffuse baryon mass Wb, 



TObh == (1 - en)mb, 



(49) 



where the energy released in electromagnetic radi- 
ation, neutrinos, kinetic energy, and possibly grav- 
itatonal radiation is 



TOomittod — <^nmh- 



(50) 



The baryon mass sequestered in massive black 
holes, TOb = rTibh/(l — £«), could be substantial 
if e„ were close to unity. However, if an appre- 
ciable part of the binding energy were released as 
electromagnetic radiation then the bounds on the 
radiation background (in category 7) and those 
of the CMBR distortion (eq. [15]) would require 
that the energy was released at very high redshift, 
kT > 10^ K. The estimates discussed in §2.7.6 in- 
dicate the mass fraction released in gravitational 
radiation is small. Thus it seems likely that e„ is 
small, as is assumed in entry 3.13. The efficiency 
factor eg typical of stellar mass black holes does 
not appear in entry 3.7, because this estimate is 
based on an analysis of the progenitor star masses. 
The estimate of the mean mass density in massive 
black holes that is used for entry 3.13 is discussed 
in §2.5.3. 

2.3.7. Intergalactic Plasma 

Entry 3.1, for the baryon mass outside galaxies 
and clusters of galaxies, is the difference between 
our adopted value of the baryon density parame- 
ter (eq. [17]) and the sum of all the other entries 
in category 3. Within standard pictures of struc- 
ture formation this component could not be in a 
compact form such as planets, but rather must 
be a plasma, diffuse enough to be ionized by the 
intergalactic radiation or else shocked to a tem- 
perature high enough for collisional ionization, but 
not dense and hot enough to be a detectable X-ray 
source outside clusters and hot groups of galaxies. 

In our baryon budget 90% of the baryons are in 
this intergalactic plasma. It is observed in several 
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states. Quasar absorption lines show matter in 
low and high atomic ionization states in the halos 
of L ~ L, galaxies, extending to radii ^ 200 kpc 
(Chen et al. 2001 and references therein). Ab- 
sorption lines of HI and Mgll reveal low surface 
density photoionized plasma at kinetic tempera- 
ture T ^ 10* K, which can be well away from 
L ^ L* galaxies, as discussed by Churchill, Vogt, 
& Charlton (2003), and Penton, Stocke & ShuU 
(2004). The latter authors estimate that 30% of 
the baryons are in this state. Absorption lines 
of O VI around galaxies and groups of galaxies 
(Tripp, Savage & Jenkins 2000; Sembach et al. 
2003; ShuU, Tumlinson, & Giroux 2003; Richter 
et al. 2003), reveal matter that may be excited 
by photoionization by the X ray background ra- 
diation and by collisions in plasma at the kinetic 
temperature T ~ 10^ K characteristic of the mo- 
tion of matter around galaxies (Cen et al. 2001). 
The detection of O VII and O VIII absorption lines 
also indicates the presence of higher temperature 
regions in the local IGM (Fang et al. 2002; Fang, 
Sembach & Canizares 2003). Improvements in the 
constraints on the amount of matter in these vari- 
ous states of intergalactic baryons will be followed 
with interest. 

For the inventory we adopt the measure of 
the concentration of dark matter around galaxies 
in equation (13), and the argument discussed in 
§ 2.3.8 that the baryons are distributed like the 
dark matter on scales comparable to the virial 
radii of galaxies. The resulting division into the 
mass in baryons near the virial radii of normal 
galaxies outside clusters (entry 3.1a), and the mass 
well away from galaxies and compact groups and 
clusters of galaxies (entry 3.1b), is presented as 
sub-components, because the sum is much better 
constrained than the individual values. 

2.3.8. Baryon Cooling 

We comment here on a simple picture for the 
cooling and settling of baryons onto galaxies. The 
sum of the baryon mass densities belonging to 
galaxies, in entries 3.3 to 3.13, is f^b.g = 0.0035. 
This is 8% of the total baryon mass. Suppose 
r^b.g consists of all baryons gathered from radius 
fg around L ^ L* galaxies, and suppose we can 
neglect the addition of baryons by settling from 
further out and the loss by galactic winds. That 
is, we are supposing that at r > r^ the ratio of the 



baryon density to the dark matter density is the 
cosmic mean value, and that the baryons closer 
in have collapsed onto the galaxies. In this pic- 
ture the characteristic radius of assembly of the 
baryons satisfies 



2ngrg(7^ 



nu 



Gp„ 



n 



b, total 



(51) 



in the limiting isothermal sphere approximation 
(eq. [9]). 

In equation (51) Ug is a measure of the number 
density of luminous galaxies. We record here our 
choices for this quantity and related parameters 
that are used elsewhere. We take 



Cr/L; ^ O.OUh^ Mpc" 



(52) 



L*B 


= 1.07x IQ^^h-^LQ, 


K-' 


-- 1.45 X lO^O/^-^Lo, 


K-' 


= 2.37 X lO^O/i^^L©, 



where Lr is the luminosity density (eq. [14]). The 
characteristic galaxy luminosity. 



(53) 



is the luminosity parameter in the Schechtcr func- 
tion, with the power law index olb — ^1-1, ctr ~ 
— 1.13, and a^ = —1.14. We refer some estimates 
of energy densities to what is known about the 
Milky Way, for which we need the effective num- 
ber density of Milky Way galaxies. In the B band 
the Milky Way luminosity is Lmw 
the effective number density is 



1.3L|j, and 



HMw = Cb/LI^w = 0.013/1^ (Mpc) 



(54) 



Almost the same density follows when referred to 
the r-band. The Milky Way parameters are dis- 
cussed in the Appendix. 

With the characteristic velocity dispersion in 
equation (10) and the characteristic galaxy num- 
ber density in equation (52) the baryon accretion 
radius defined by equation (51) is 



30h-^ kpc. 



(55) 



at density contrast 1.3 x 10"* and plasma den- 
sity Ugasirg) ~ O.OOT/i^ cm~^. If plasma at this 
radius were supported by pressure at the one- 
dimensional velocity dispersion a ~ 160 km s~^ 
in equation (10) the temperature would be T ~ 
2 X 10® K. At this density and temperature the 
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thermal bremsstrahlung cooling time would be 
short enough, ~ 4 x 10^ yr, that stars would have 
formed and disks matured at z ^^ 1. 

We have lower bounds on the cooling radius 
from the observation that the neutral atomic hy- 
drogen density around L* galaxies reaches A^hi = 
1.8 X 10^" cm^^ at the effective radius ~ 20h^^ kpc 
(Bosma 1981), and Mg II absorption lines are ob- 
served at radius r^ 40ft.^^ kpc (Steidcl, Dickinson 
& Persson 1994). For the Milky Way the distri- 
bution of RR Lyr stars cuts off sharply at 50 kpc 
(Ivezic et al. 2000). The cooling radius must be 
larger than these indicators of relatively cool mat- 
ter. Equation (55) is not inconsistent with this 
condition. Though the history of baryon accre- 
tion by galaxies undoubtedly is complex, we can 
imagine, as a first approximation, that a substan- 
tial fraction of the baryons now concentrated in 
galaxies are there because they were able to cool 
and settle from an initial distribution similar to 
that of the dark matter. 

The relative distributions of baryons and dark 
matter at distances much larger than Vg from 
galaxies might not be greatly disturbed from the 
primeval condition. If so, then the product of the 
baryon density parameter with the virialized dark 
matter mass fraction in equation (13) is an esti- 
mate of the baryon mass that resides within the 
virial radii of normal galaxies, and the remainder. 



n 



b,ig 



0.016, 



(56) 



which is presented as entry 3.1b, would be located 
outside galaxies and remain less than fully docu- 
mented. 

2.4. Primeval Gravitational Energy 

In the ACDM cosmology the gravitational bind- 
ing energy of the present mass distribution has two 
contributions. The first, which we term primeval, 
is a result of the purely gravitational growth of 
mass fluctuations out of the small adiabatic de- 
partures from a homogeneous mass distribution 
present in the initial conditions for the Friedmann- 
Lemaitre cosmology. The second, to be discussed 
in the next subsection, is the result of dissipa- 
tive settling of baryons that produced the baryon- 
dominated luminous parts of the galaxies along 
with stars and star remnants. We can find sensi- 
ble approximations to the primeval and dissipative 



components because, as we will discuss, the char- 
acteristic length scales are well separated. 

The primeval gravitational energy is defined by 
imagining a universe with initial conditions identi- 
cal to ours in all respects except that the baryonic 
matter in our universe is replaced by an equal mass 
of CDM in the reference model. At the present 
world time this reference model contains a clus- 
tered distribution of massive halos with gravita- 
tional binding energy density that we term the 
primeval component. ^^ This component is esti- 
mated as follows. 

The Layzer (1963) - Irvine (1961) equation for 
the evolution of the kinetic and gravitational ener- 
gies of nonrclativistic matter, such as CDM, that 
interacts only by gravity is 

— {K + W) + -{2K + W)=0. (57) 
at a 

The kinetic energy per unit mass is 

K = {mv^/2)/{m), (58) 

where v is the peculiar velocity of a particle with 
mass m. The gravitational potential energy per 
unit mass is 



W = -^Gp„ 



cPr£,{r)/r, 



(59) 



where pm is the mean mass density and ^ is the 
reduced mass autocorrelation function. 

We can use a simple limiting case of the Layzer- 
Irvine equation, because in the ACDM cosmology 



^^Onc surely would say that in this model universe the viri- 
alized dark matter halos have gravitational binding energy. 
Since there was no energy transfer to some other form, one 
might also want to say that this binding energy must have 
been present in the initial conditions. Furthermore, one can 
assign to a linear mass density fluctuation with contrast 
5{t) > 0, and comoving radius x (physical radius xa(t)), 
a gravitational energy per unit mass, W' ~ —G{p)5{ax)'^, 
which is constant in linear perturbation theory and com- 
parable to the binding energy of the final virialized halo. 
Perhaps one can use this as a guide to a definition of the 
primeval energy belonging to the density fluctuation, de- 
spite the problem that the mean of W vanishes, and the 
fact that in general relativity theory there is no general 
definition of the global energy density of a statistically ho- 
mogeneous system. We have not been able to find a useful 
approach along these lines. We might add that if our uni- 
verse had been Einstein-de Sitter then we would have de- 
fined the primeval energy density by a numerical solution 
of equation (57). 
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the universe has now entered A-dominated expan- 
sion, which has caused a significant suppression of 
the rate of growth of large-scale departures from 
homogeneity. This means that the first term in 
equation (57) has become small compared to the 
second term. Thus it is a reasonable approxima- 
tion to take 2K = —W, the usual virial equilib- 
rium relation. Then the gravitational binding en- 
ergy per unit mass is U — K + W — W/2, and 
the cosmic mean primeval gravitational binding 
energy is the product of U with the mean mass 
density fi-m in matter (eq. [8]). With the normal- 
ization P{k) — j (Pr^{r)e^^''^' for the mass fluc- 
tuation power spectrum, the expression for the 
primeval gravitational binding energy in this ap- 
proximation is 



n 



BE,p 



167r2 



dkP{k). (60) 



For a numerical value we use the mass fluctua- 
tion power spectrum P(k) in Figure 37 of Tegmark 
et al. (2004b) at fc < 0.1ft, Mpc"^ At smaller 
scales the Tegmark et al. spectrum decreases too 
rapidly with increasing wavenumber fc to be a good 
approximation to the present mass distribution 
(Davis & Peebles 1983). We approximate the spec- 
trum as 



P{k) = 7000 



O.l/i Mpc" 



1.23 



h-"^ Mpc^ (61) 



at k > O.lh Mpc^^. This is based on the Fourier 
transform of the pure power law model for the 
galaxy autocorrelation function (eq. [5]), 

P = 4:TTr2k'''-^r{2 - 7) sin 7r(2 - 7)/2, (62) 

which gives P = 5200/i~^ Mpc^ at fc = O.l/i Mpc"^ 
We choose the somewhat larger normalization in 
equation (61) to fit the SDSS measurement. The 
numerical result is 



P{k) dk = 4800h-^ Mpc^ 



(63) 



The integral of the Tegmark et al. (2004b) power 
spectrum over all wavenumbers is 2700/i~2Mpc, 
which is 0.55 times the value in equation (63). Our 
estimate of the integral is based on measurements 
of the actual power spectrum, not the spectrum 
that would have obtained if there were no dissi- 
pative settling of baryons, but the error is small 



because the integral is dominated by lengths large 
compared to the scale of separation of baryons 
from the dark matter. 

Equations (60) and (63) yield the entry in line 
4 of Table 1, 



n„ 



-7.7 X 10-"^. 



(64) 



This is the density parameter of the gravitational 
binding energy of the present departure from a 
homogeneous mass distribution, ignoring the ef- 
fects of the dissipative settling of baryons. It will 
be useful to note that a measure of the length 
scale of the gravitational energy is the half-point 
of the integral in equation (63), at wavenumber 
ki/2 = 0.27 h Mpc^^, or half wavelength 



Al/2 = 7r/fci/2 



Uh^^ Mpc. 



(65) 



The kinetic energy per unit mass belonging to 
equation (64) is 



K 



-[/ = 3(410 km s"')72 



(66) 



The velocity in parentheses is the one-dimensional 
single-particle line-of sight rms peculiar velocity. 

The gravitational binding energy is not equal to 
a sum over the contributions from individual ob- 
jects, but we can write useful approximations to 
the decomposition into the three components — 
the virialized parts of the massive halos of L ~ L* 
galaxies, the rich clusters, and large-scale cluster- 
ing — shown in category 4 in the inventory. 

Since galaxy rotation curves tend to be close to 
flat we write the binding energy of the virialized 
parts of the dark matter halos of the galaxies as 



n. 



Pn 



10 



-7.2 



(67) 



where ilm is the matter density parameter (eq. [8]), 
a is the characteristic velocity dispersion in equa- 
tion (10) and the last factor is the virialized mass 
fraction (eq. 13). This is about 10% of the to- 
tal gravitational binding energy. The dissipative 
settling that produced the baryon-dominated lu- 
minous parts of the galaxies would have perturbed 
the massive halos, but the disturbance to the 
primeval gravitational energy is small because the 
baryon mass fraction is small. 

Our estimate of the binding energy of the dark 
halos of the galaxies may be compared to the den- 
sity parameter given by equation (60) when the 



16 



integral over P{k) is restricted to small scales, 
k > Tr/ry, where r^ is the virial radius in in equa- 
tions (10) and (67). The result, with equation (61) 
for P(fc), is twice the value in equation (67). The 
difference is an indication of the ambiguity of sep- 
arating gravitational energy into components. For 
entry 4.1 we adopt the value in equation (67) as 
the more directly interpretable. 

Equation (67) is a reasonable approximation for 
most of the mass in galaxy-size dark halos in lu- 
minous field galaxies such as the Milky Way, but 
in rich clusters the galaxies tend to share a dark 
halo that is close to smoothly distributed across 
the cluster. Our estimate in entry 4.2 for the 
primeval gravitational binding energy belonging to 
rich clusters follows equation (67), with a — 800 
km s""'^ and fl^i from equation (46). Rich clusters 
share about 15% of the total gravitational binding 
energy. 

Entry 4.3 is the difference between the total in 
entry 4 and the sum of entries 4.1 and 4.2. The 
Layzer-Irvine equation indicates that the binding 
energy is dominated by a length scale (eq. [65]) 
that is much larger than galaxy virial radii. Con- 
sistent with this, entry 4.3 is larger than entry 4.1 
for the binding energy of the dark halos of galax- 
ies. The difference is not large, however, because 
at small scales the integral over the power spec- 
trum converges slowly, as k~'^'^^. 

2.5. Dissipative Gravitational Settling 

Dissipative settling has increased the magni- 
tude of the gravitational binding energy from that 
prescribed by the primeval conditions considered 
in the last section. In §2.5.1 wc discuss the energy 
released in producing the increased mean density 
of baryons relative to dark matter in the lumi- 
nous parts of the galaxies, in §2.5.2 wc estimate 
the gravitational energy released in stellar forma- 
tion and evolution, and in §2.5.3 we consider the 
central massive compact objects in galaxies. 

2.5.1. The Luminous Parts of Galaxies 

In the Milky Way galaxy the mass within our 
position, at about 8 kpc from the center, is roughly 
equal parts baryonic and dark matter, or about 6 
times the cosmic mean ratio (cq. [19]). This is 
thought to be the usual situation in the luminous 
parts of normal galaxies. The amount of gravi- 



tational binding energy released in producing this 
concentration of baryons depends on how it was 
done. In one limiting case one may imagine that 
stars formed in the centers of low mass dark halos 
with relatively small dissipation of energy — apart 
from star formation — because the depths of the 
gravitational potential wells were small, and that 
the low mass halos later merged without any ad- 
ditional dissipation, the dense baryon-dominatcd 
parts remaining near the densest regions to form 
the present-day baryon-dominated luminous parts 
of galaxies. (This is an extreme version of the 
scenario discussed by Gao et al. 2003). In an- 
other extreme, one may imagine that the baryons 
settled into previously assembled galaxy-scale ha- 
los, which would dissipate considerably more en- 
ergy. A galaxy has a definite computable gravi- 
tational binding energy, of course (apart from the 
difficulty of correcting for ongoing accretion), but 
to relate this to the energy dissipated in producing 
the galaxy would require an analysis of what the 
mass distribution would have been in the absence 
of dissipation, which is not an easy task. 

These considerations lead us to offer only a 
crude estimate for entry 5.1, which we write as 
the product of the density parameter belonging to 
baryons in galaxies — the sum rib,g — 0.0035 of 
the density parameters in entries 3.3 to 3.13 — 
with the kinetic energy per unit mass, K — 3(T^/2 
and a = 160 km s^^. The result is a 2% ad- 
dition to the primeval halo gravitational binding 
energy (entry 4.1). If the baryon concentrations 
in galaxies formed at high redshifts in small halos 
the dissipative energy released would be an even 
smaller fraction of the total. 

2.5.2. Stellar Binding Energy 

The amount of binding energy released in star 
formation is easy to define because the relative 
length scale is small. We write the gravitational 
binding energy per unit mass for a star with mass 
m and radius r as 



BE 

m 



-K 



Gm 



(68) 
1.74, and K ^ 



The prefactor for the Sun is Kq 
0.3 for a homogeneous sphere. 

For main sequence stars we use the zero-age 
mass-radius relation, r ~ 0.85 to"-^" for 0.08 < 
TO < 0.79, r 2± 0.93 to!" for 0.79 < m < 1.38, 
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and r ~ 1.15 m"'^^ for 1.38 < m < 100, in so- 
lar units. These numbers are assembled from Ezer 
& Cameron (1967), Cox & GiuUi (1968) and Cox 
(2000). Integration of Gm^/r over the PDMF 
gives BE/m = 3.7 x 10"*^. The product of the last 
number with the density parameter of the mass 
in main sequence stars (entries 3.3 plus 3.4), with 
K c:i Kq , is the estimate of the gravitational bind- 
ing energy, JIbe = — lO"^'"'^, for stars. We similarly 
obtain the substellar gravitational binding energy, 
f^BE = — 10~^-^, where r is fixed at 0.096 r© (Bur- 
rows et al. 2001). This is a small addition to the 
sum in entry 5.2. 

We construct a model for the white dwarf mass 
function from an approximation to the relation be- 
tween the progenitor main sequence mass and the 
white dwarf remnant mass (Claver et al. 2001; 
Weidemann 2000), 



TOwd ~ O.OSmnis + 0.45 m0. 



(69) 



and our IMF. White dwarf masses run from 
0.53 to 1.09 ?Ti0 for the main sequence mass 
range 1 < m^g < 8 ttiq we have adopted. 
The white dwarf mass function dN/dnivjd = 
{dN/dmais)id'm^s/d'm^d) thus obtained agrees 
well with the observed white dwarf mass distri- 
bution of Bergeron & Holberg for m > 0.5™© (in 
preparation). Here we ignore low mass helium 
core white dwarfs. In our mass function the mean 
white dwarf mass is 



("Iwd) 



0.62 m0. 



(70) 



From our mass function and the mass-radius re- 
lation given by Shapiro and Teukolsky (1983) we 
obtain the mean white dwarf gravitational binding 
energy per unit mass. 



BE 

m 



-58K. 



wd~ 



Gm(7 



-1.2 X 10" 



(71) 



Since the fractional half-mass radius of a white 
dwarf is 0.57 times the solar value (Schwarzschild 
1958), we have taken Ky^^ ~ 1.0. The product 
with the mass density in white dwarfs (entry 3.5) 
is entry 5.3. 

We take the binding energy of a neutron star 
to be 3 X 10^"^ erg (e.g., Burrows 1990; Janka & 
Hillebrandt 1989), or 



The product with entry 3.6 is entry 5.4, Obe,ns — 
10~^'^. The gravitational binding energies in neu- 
tron stars and stellar mass black holes are sub- 
stantially larger than the gravitational binding en- 
ergies in all other forms. 

2.5.3. Black Hole Binding Energy 

Our definition of the binding energy associated 
with a black hole requires careful explanation be- 
cause it has some curious properties, including vi- 
olation of the thought that it would be logical to 
consider the mass of a black hole to be purely grav- 
itational if the matter out of which it formed has 
lost its existence. 

We choose the definition by analogy to nuclear 
and Newtonian gravitational binding energy, in 
terms of the energy liberated in the assembly of 
a system out of its initial parts, that is, the dif- 
ference between the total mass of the initial parts 
and the mass of the assembled system. In the same 
way, we use equations (49) and (50) to define the 
binding energy of a black hole by the difference be- 
tween the mass rrib of the initial parts — baryons 
— and the mass Wbh = (1 — in)'iTih of the final 
black hole. Thus our definition of the binding en- 
ergy of a black hole is 



BE = -e„mb = - 



1 



-mbh- 



(73) 



BE/m 



-0.12. 



(72) 



The magnitude of BE is the energy emitted as 
electromagnetic and gravitational radiation, neu- 
trinos, and kinetic energy, as is appropriate for our 
purpose of telling the energy transfers and balanc- 
ing the baryon budget. It will be noted that in 
this definition the binding energy depends on how 
the black hole formed. For example, a Solar mass 
black hole that formed with efficiency e ~ 0.99 
is assigned binding energy —99mQ, because it re- 
leased that much energy, while an identical black 
hole that formed with e„ — 0.01 is assigned a very 
different binding energy, — O.OItti©. 

Entry 5.5 for the gravitational binding energy 
of stellar mass black holes is the product of entry 
3.7, which is our estimate of the baryonic mass 
entering the black hole, with the efficiency fac- 
tor £„. In the standard picture for the formation 
of a stellar mass black hole, a core of baryons is 
first burned to heavy elements, and the subsequent 
collapse to a black hole may release little more en- 
ergy. In this case the efficiency factor could be as 
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small as £„ ^ 0.009, which is the binding energy 
released as starlight. It could also be as large as 
Cs ~ 0.03 if the collapse proceeded through a pro- 
toneutron star as an intermediate state. It can- 
not be much larger, however, without violating 
the constraints from the radiation energy density 
(see §2.7) and the relic supernova neutrino flux at 
Super-Kamiokande (Fukugita & Kawasaki 2003). 
One way to estimate the mass density in the 
massive black holes in the nuclei of galaxies uses 
the correlation of the black hole mass with the 
bulge luminosity. A convenient approximation to 
the relation, for B-band luminosities, is (Gebhardt 
et al. 2000; Ferrarese 2002; see also Kormcndy & 
Richstone 1995) 



The numerical result, 



M,/m(. 



10 



-2.0±0.3 



L 



bulge/ -tJ© 



/Lq, 



(74) 



FHP estimate that the fraction of the B-band lu- 
minosity density in ellipticals and SO galaxies is 
0.24, and the fraction in the bulges of spheroids 
is 0.14. The products of equation (74) with the 
luminosity fractions and the luminosity density in 
equation (14) gives the mass density parameters 
in massive black holes. 



f}. (early) = 10 
f}.(late) = 10 



-5.6±0.3 



-5.9±0.3 



(75) 



Salucci et al. (1999) give a consistent, but slightly 
larger value. 

For early-type galaxies we can use the tight re- 
lation between the black hole mass and the bulge 
or spheroid velocity dispersion (Merritt & Fer- 
rarese 2001; Tremaine et al. 2002). The Sheth 
et al. (2003) estimate of the velocity dispersion 
function for early-type galaxies is 



dN ^4>^ 



^-e-(^/-)^ (76) 



with a = 6.5, (3 — 1.93, Gi, — 89 km s ^, 
(/), = 0.0020 Mpc-^ The Tremaine et al. (2002) 
estimate of the black hole mass- velocity dispersion 
relation is 

M. = B{<j/<JhT. (77) 

with B = 1.3 X lO^m©, a = 4.0, and cr^ = 200 
km s^^. The product of the two expressions, inte- 
grated over CT, gives the mean mass density. 



T{{a + a)/l3) / a^ 
' T{a/(3) [ah 



(78) 



17, (early) = 10" 



(79) 



is close to but smaller than the more direct esti- 
mate in equation (75). Although the formal un- 
certainty in equation (79) is smaller it rests on the 
condition that the Sheth et al. galaxies are a fair 
sample of the early-type galaxies, which will re- 
quire careful debate. "'^^ Thus in the inventory we 
quote equation (75). 

2.5.4- Quasar Luminosities and Remnants 

Soltan (1982) and Chokshi & Turner (1992) 
have considered the relation between the rate of 
radiation of energy by quasars and AGNs and the 
accumulation of mass in the quasar engines, which 
are assumed to be massive black holes in the cen- 
ters of galaxies. In this repetition of the calcula- 
tion we take the number of quasars per unit lumi- 
nosity and comoving volume to be 



L, 



L,$(i, 



dn 



dL {L/L^)^ + {L/L^)f^' 



(80) 



where, from Groom et al. (2004), a = 3.31, (3 = 
1.09, 

L,$(L,) = 1.81 X 10"^ Mpc-^ (81) 

and the present characteristic luminosity is 

i, =6.7x lO^°iij(0). (82) 

In the Groom et al. luminosity evolution model 
this luminosity evolves as L^,{z) ex iQ^-39z~a.29z 
to redshift z = 2.1, the deepest rcdshift used in the 
Groom et al. analysis. The peak of the observed 
comoving quasar number density is at z ~ 2.5, 
and at 2.5 < 2; < 5 the density varies about as 
n (X e-i-5^ (Fan et al. 2001, Fig. 3). As a con- 
venient approximation to this behavior we adopt 
the Groom et al. (2004) luminosity evolution at 
z < 2.1, constant comoving luminosity density 
from z = 2.1toz = 3, and negligibly small lu- 
minosity at larger redshifts. 

The integral J dL L dn/ dL is the comoving lu- 
minosity density, and the time integral multiplied 



l^The velocity function of Sheth ct al. gives (cr"')!/-' = 180 
km s~^, compared to our estimate of the characteristic ve- 
locity dispersion, a-» = 200 — 220 km s~^, in early-type 
galaxies. Perhaps this is related to the difference. 
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by the bolometric correction is the net comov- 
ing density of energy released. Using the Elvis 
et al. (1994, Table 17) bolometric correction fac- 
tor BC = Lx,o\l{vL^)\^^^^^ = 12, and ignoring 
the difference between B and bj passbands (since 
the quantity that concerns us, vL^, for quasars, 
is close to flat), we estimate that the integrated 
energy density released by the quasars is 



^ 



QSOcm 



= 8 X 10" 



This uses the solar luminosity, 

vL^{Q) = 2.22 X 10^^ erg s 



(83) 



(84) 



at A = 4450 A. 

Before comparing this estimate to the accu- 
mulated mass in black holes let us check consis- 
tency with the integrated background radiation. 
Equation (83) with the Elvis et al. bolomet- 
ric corrections indicates that the integrated back- 
ground from quasars at 100/i < A < 1000 A is 
ri ~ 3 X 10"*, or 1% of the total (eqs. [106] 
plus [107] below). To estimate the expected X- 
ray background we imagine the radiation from the 
quasars all is emitted at effective redshift z = 2. 
In this simple model the present energy density 
per logarithmic interval of frequency at 2 keV is 



vVt^{2 keV) = 



^ 



QSOcm 



VpL^ 



1 



+ z J^ L^dv ' 



(85) 



where on the right-hand side hv^. — 2(1 + z) keV. 
The mean spectrum in figure 10 in Elvis et al. 
(1994) for radio-quiet quasars, with the bolometric 
factor, indicates that v^L^^ = J^ L^dv/bQ at rest 
frame energy hv^ = 6 keV. These numbers give the 
present energy density per logarithmic interval of 
photon energy vVli, = lO^^"^ at 2 keV. The mea- 
sured value of the X-ray background is (De Luca & 
Molendi 2004) j/fJ^ = lO'^-* at 2 keV, about three 
times what is indicated by equations (83) and (85). 
Since about 80% of the X-ray background at 2 keV 
is resolved (Mushotzky et al. 2000; Worsley et al. 
2004), these numbers allow room for a significant 
population of optically faint quasars. 

We turn now to the efficiency e„ for production 
of electromagnetic radiation in the accumulation 
of the present mass in quasar remnants. If e„ is 
small the estimate of the integrated mass added 
to the black holes by the observed energy produc- 



(eq. [83]). The ratio of this expression to the mass 
density in massive black holes (the sum of entries 
5.6 and 5.7) is our estimate of the radiation effi- 
ciency, 

e„ = 0.02. (86) 

This is one fifth of the commonly discussed value, 
e„ ~ 0.1. Since our estimate of the X-ray back- 
ground from optically identified quasars is one 
third of the measured value, it may be that equa- 
tion (83) is low by a factor of about three and 
equation (86) accordingly low by a like factor. A 
closer check of consistency of the idea that the 
massive black holes in the centers of galaxies are 
the quasar remnants awaits advances in surveys of 
optically faint quasars in broader ranges of wave- 
length and redshift and better understanding of 
the quasar emission mechanism. We may also 
hope that future work will establish the natures of 
the sources of the harder X-ray background and 
the possible relevance to the black hole mass bud- 
get. 

2.6. Nuclear Binding Energy 

2.6.1. Heavy Element Abundances 

We consider here the binding energy released by 
nuclear burning in stars. We normalize the heavy 
element abundances to the Solar mass fractions in 
hydrogen, helium, and heavy elements, 

^0 = 0.71, Yq = 0.27, Zq ^ 0.019. (87) 

The ratio Z/X =- 0.027 is derived (BP2000) as the 
initial solar value from (Z/X)soiar surface — 0.0230 
of Grevesse & Sauval (2000), and X is also the ini- 
tial value of BP2000. The metallicity in star pop- 
ulations is correlated with the galaxy luminosity. 
An average over the Schechter luminosity function 
of the Kobulnicky & Zaritsky (1999) correlation of 
the oxygen abundance with the B-band galaxy lu- 
minosity, taking 12-f'log(0/H)Q = 8.83 (Grevesse 
& Sauval 2000)^'^ as the zero point, indicates that 



tion by quasars and AGNs is Af^BH — ^ 



QSOcm 



h-n 



^^Recent work on the solar heavy element abundance 
suggests a significantly lower oxygen abundance, 12 -|- 
log(O/H)0 = 8.69 (AUende Prieto, Lambert & Asplund 
2001), but if the heavy element abundances of other ele- 
ments are scaled down in a similar manner, as indicated 
by the same team, our net result is not affected. In fact, 
Bahcall and Pinsonncault (2004) give {Z j X^^aXm surface = 
0.0176 based on the new abundance, which leaves our result 
unchanged. 
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Table 3 
Heavy Element Masses 



Objects 


Mean Metallicity 


Composition 


lO^Oz 


1 main scq. stars 


(Z) 


solar 


3.2 


2 substellar objects 


(Z) 


solar 


0.2 


3 white dwarfs 


1 


c+o 


36 


4 cool gas 


{Z) 


solar 


1.3 


5 clusters 


Z0/3 


solar 


1.1 


6 warm plasma 


Z0/3O 


solar 


2.5 


7 neutron stars^ 


1 


(Fe) 


5 


8 stellar mass black holcs^ 


1 


(Fc) 


6.8 


9 massive black holes^ 


Ze 


(solar) 


0.005 


^vanished 









the mean metallicity in galaxies is 

{Z) = 0.83(1 ± O.3)Z0 = 0.016 ± 0.003. 

Table 3 lists our estimates of the density param- 
eters belonging to the mass in heavy elements in 
several categories of objects. The entry for main 
sequence stars in the first line is the product of {Z) 
with the sum of the density parameters in entries 
3.3 and 3.4 in the inventory, and the second line 
uses entry 3.8. Since the main elements in white 
dwarfs are carbon and oxygen, with thin hydro- 
gen and/or helium layers that typically amount to 
< 0.1% of the mass, we enter in the third line the 
density parameter from entry 3.5. We assign the 
cool gas in entries 3.9 and 3.10 the same mean 
metallicity as the stars (eq. [88]). The metallic- 
ity of the intracluster plasma (entry 3.2) is ob- 
served to be about one third of solar (Mushotzky 
& Loewenstein 1997; Fukazawa et al. 1998; White 
2000), as indicated in Table 3. We suppose the in- 
tergalactic plasma (entry 3.1) may have metallicty 
about 3% of Solar. This likely is larger than the 
metallicity in plasma in the voids (Penton, Stocke 
& ShuU 2004), and smaller than the metallicity 
in the plasma observed as absorption line systems 
around galaxies (Sembach et al. 2003; Churchill, 
Vogt, & Charlton 2003), and perhaps is a reason- 
able factor-of-three compromise. We assume that 
neutron stars and stellar mass black hole form by 
the collapse of an iron core, and that massive black 
holes grew by the accretion of matter with about 
the Solar heavy element abundance. These heavy 
elements are entered in Table 3, but they are now 
sequestered from the inventory. 

Our estimate of the total production of heavy 
elements, including those that have been lost in 



neutron stars and black holes, is 

Viz = 5.7 ±1.2 X lO"''. 



(89) 



The heavy elements in white dwarfs amount to 
about 65% of the total. The matter in this large 
reservoir is liberated only on the rare occasions of 
Type la supernovae. 

The model for the rate of Type la supernovae is 
uncertain; we consider the widely accepted Whe- 
lan & Iben (1973) binary white dwarf picture and 
use a simple model for the supernova rate. 



i?SN Ia(t) = A / dtV(i') 



8mo ^^ 

dm— — X 

.j^.^ dm 

exp(-[t-f'- At]/T-), (90) 

where A is the normalization determined by the 
empirical Type la supernova rate at zero redshift, 
which is obtained from the three surveys men- 
tioned earlier (eq. [38]), 

i?SNia(to) = 0.0027+°°°!^ (100yr)-iMpc-3, 

(91) 
and At{mi,) — 13(m/TO0) ^'^ Gyr -I- 5 is the time 
for the formation of white dwarfs plus the time 
delay to form the Roche-lobe contact, and we take 
(5 = to 1 Gyr. The minimum mass is chosen to 
be mniin = max[3m0, {(t — t')/13Gyr}^"''], where 
3m0 corresponds to the white dwarf mass O.7to0 
according to eq. (69) so that the Chandrasekhar 
limit is met. This is the model taken by Madau 
et al. (1998) (see also Gal- Yam & Maoz 2004). 
We adopt t — A Gyr to account for a moderate 
increase (by a factor of 2 ± 1) of the observed SNIa 
occurence to z = 0.6 (Pain et al. 2002). i'' The 



'A short time scale is in conflict with the Type la super- 
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effective time span for the cumulative occurrence 
of Type la supernovae normalised to the rate at 
2; = is 



Tcs 



-fi'sNIa(^o) 

15 ± 3 Gyr, 



to 



RsNIa{t)dt 



(92) 



for our model of star formation (eqs. [31] and [33] 
with Zf = 2.5 to 4). The cumulative comoving 
number density of supernova in this model is 4.0 x 
10^ Mpc~^. Since we estimate that the number 
density of white dwarfs created is 8 x 10^ Mpc""^, 
this would mean that about 1% of the white dwarfs 
have been disrupted. If OJitiq of '"'^Ni is produced 
in each supernova (e.g., Branch 1992), the mass 
density of iron-group elements produced by type 
la supernovae is rjpo = 2.0 x 10^®. This is supple- 
mented with Ope — 3.6 x 10^^ from Type Il+Ib/c 
supernovae (see eq. [39]), which produces O.OTStti© 
of iron per event (Arnett 1996; Weaver, Zimmer- 
man & Woosley 1978). These estimates indicate 
that the total density parameter in the iron group 
elements is 

r^Fc = 6 X 10"^ (93) 

of which about 60% is from Type II supernovae. 

This estimate can be compared to the product 
of the mass in heavy elements not locked up in 
stellar remnants (the sum of lines 1,2, and 4 to 6 
in Table 3), 



n'r, ^ 0.8 ± 0.25 X 10" 



(94) 



with the iron group mass fraction (Grevesse & 
Sauval 2000), 



f^Fc = Omin'r, ^ 6.3 X 10" 



(95) 



The two approaches give a consistent picture for 
the origin of iron and the supernova rates. 



nova rates in early type galaxies. With t = 4 Gyr, the 
SN fraction in early type galaxies is 0.30 at z Si 0, which 
is consistent with the observed value, 0.35_g'j^g. Here we 
have identified stars formed at 2: > 0.7 as an early popula- 
tion that is partitioned into early-type galaxies and bulges 
of disk galaxies according to the FHP fractions of bulge lu- 
minosities. The fraction drops to < 0.25 if r = 3 Gyr. The 
observed rates of Type la supernovae in morphologically 
separated galaxies seem to be proportional to the r band 
luminosity density; the luminosity density in E/SO galaxies 
is 31% of the total (Nakamura et al. 2003). This would 
suggest a close to constant supernova rate. 



The growth of the abundance of heavy elements 
is accompanied by the accumulation of helium, 
apart from the heavy elements that enter white 
dwarfs, neutron stars and stellar mass black holes. 
The estimate of AY/AZ in equation (20), applied 
to all entries in Table 3 except 3, 7 and 8, is our 
estimate of the present helium mass fraction in 
excess of primeval, 

AY 
Any ^ -r— ^z = 1.7 ± 1.0 X 10"*. (96) 

This is consistent with our estimate from the prod- 
ucts of stellar evolution (eq. [120]). 

There are in the literature many analyses of 
the stellar production of heavy elements (see for 
example Calura & Matteucci 2004 and references 
therein). Our results are generally consistent for 
relevant entries. Important differences are our in- 
clusion of heavy elements in white dwarfs, which 
dominate the entry, and our consideration of iron 
in the progenitors of supernovae. These compo- 
nents are crucial to the consideration of the bal- 
ance between the nuclear binding energy stored in 
the elements and the energy released in radiation. 

2.6.2. Nuclear Binding Energies 

Our definition of the nuclear binding energies 
of the heavy elements differs from the usual prac- 
tice in tables of nuclei. Because we are interested 
in the release of energy in the formation of the 
heavy elements, we calculate the binding energy 
with respect to free protons and electrons. 

We write the energy released in the production 
of the heavy elements present in the interstellar 
medium, and in stars when they formed, as 

^NB, z = (0.0081 + 0.0071 AYlAZ)nz, (97) 

where 0.0081 is the energy generation efficiency 
factor for the solar composition. The nuclear bind- 
ing energy in substellar objects and in diffuse mat- 
ter, in entries 6.2, 6.4, and 6.5 in Table 1, is com- 
puted from this equation with the heavy element 
masses in Table 3. For main sequence stars, we 
add to equation (97) the nuclear binding energy 
associated with the helium that has been produced 
in the stars, which we take to be on average 5% 
of the star mass, that is, half the helium a star 
produces while it is on the main sequence. The 
sum is our estimate of the nuclear binding energy 
in stars in entry 6.1. 
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Entry 6.3, for white dwarfs, is 



2.7.1. The X ^ Ifi Background 



n 



NB,wd — 



[0.0080(1 - Yp) 
10"^-^ 



O.OOlOlplf^z 



(98) 



The second term in parentheses accounts for the 
primordial heKum abundance (eq. [18]). The value 
of riz in the first line is taken from line 3 of Ta- 
ble 3. This component amounts to about 40% of 
the nuclear binding energy. The large amount of 
nuclear burning in the stellar giant and supergiant 
phases is discussed in connection with Table 5 be- 
low. 

The nuclear binding energy in the matter out 
of which a neutron star formed was converted to 
gravitational binding energy by the dissociation 
of the heavy elements during a supernova. The 
nuclear binding energy released in the formation 
of the heavy elements that were part of the raw 
material for a neutron star added to the radiation 
background, of course, but that accounting now 
belongs in category 5. The total entered in the 
inventory for category 6 accordingly is about 15% 
smaller than the sum of the nuclear binding ener- 
gies in all entries in Table 3. The nuclear bind- 
ing energy to compare to the energy required to 
produce the radiation background is the full sum, 
including iron core progenitors for neutron stars 
and black holes. 



n 



NB 



-5.7 ± 1.3 X 10"^ 



(99) 



We discuss the relation to the energy density in 
radiation in the next subsection. 

2.7. The Radiation Backgrounds 

The cosmic energy density in electromagnetic 
radiation is thought to be dominated by mildly 
redshifted starlight, at A '^ 1/i, and a far infrared 
peak at A '^ lOOfi that is produced by the ab- 
sorption and reradiation of starlight and the light 
from AGNs (Hauser & Dwek 2001 and references 
therein). The energy densities at radio and X-ray 
to 7-ray wavelengths are much smaller, but they 
are useful measures of high energy processes, as is 
the gravitational wave background. Neutrino pro- 
duction might be counted as part of the radiation 
backgrounds, but we find it convenient to enter 
neutrinos in a separate category. 



Observations of the optical to near infrared 
extragalactic background light that report posi- 
tive detections are summarized in Table 4. They 
suggest that the surface brightness per logarith- 
mic interval in frequency is about constant at 
i^I^ = 20 ± 5 nW m"2 si~^ in the range 3500 A to 
3.5 /i. This corresponds to energy density 



a 



opt 



2.3 ±0.6 X 10" 



(100) 



in the optical to near infrared. In view of the 
technical difficulty of these observations, equa- 
tion (100) may conservatively be taken as an upper 
limit. Integrated galaxy number counts give sur- 
face brightness typically one third of the entries in 
Table 4 (Madau & Pozzetti 2000; see also Hauser 
& Dwek 2001, Table 3), 



a 



opt 



0.9 ±0.2 X 10" 



(101) 



This might be considered a lower limit. 

We have a check from the energy density com- 
puted as a time integral of the luminosity density. 



dX 
T 



dt 



1 



■X{t)C 



A(t) 



(102) 



The integral is from high redshift to the present 
world time, to. The integrand is the luminosity per 
logarithmic interval of wavelength and comoving 
volume, corrected by the redshift factor, I + z = 
l/a{t), and evaluated at the redshifted wavelength 
X{t) = a(t)Ao, where the observed wavelength is Aq 
and the expansion parameter a{t) is normalized to 
unity at the present epoch. In the flat cosmological 
model the integral is 



u ^ H^ ^XoCxJ, 



/ = 



da 



y/n^m/a^ + i-nxJ KCx 



dXC{X,z) 



(103) 



We take the shape of the present cosmic spectrum 
XL\ as a function of wavelength from Figure 13 
in Blanton et al. (2003), and we normalize to the 
luminosity density at A '--^ 1/z in equation (14). We 
approximate the evolution of the comoving lumi- 
nosity density by extrapolating the Rudnick et al. 
(2003) power law fits to the evolution in the rest 
frame U, V and B bands, in the form 

/:(A,z)cx(l + z)'^, /3 = 0.93(4400A/A)2-5. (104) 
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Table 4 
Radiation Background 



A (m) 


ul„ (nW m"^ sr"^) 


0.10 


0.60" 


0.30 


13 is*" 


0.37 


24 ±8" 


0.55 


18 ±8'" 


0.81 


24 iQ*" 


2.2 


23 ig"" 


3.5 


12 ±3"^ 



'Henry (1999) ''Bernstein ot al. (2002) 
"Matilla (1990) ''Wright & Rccsc (2000) 



We truncate the integral at z = 3, the hmit of 
the Rudnick et al. measurements. In this model 
the dimensionless integral in equation (103) is / = 
1.35. (If the comoving luminosity density were 
constant the integral would be / = 0.82.) The 
result is 



r!opt = 1.5 X 10-^ 
For the inventory we adopt 



f^opt = 1.6 ± 0.7 X 10" 



(105) 



(106) 



The error spans the estimates based on mea- 
surements of the surface brightness of the sky 
(eqs. [100]), the source counts (eq. [101]) and the 
luminosity density (eq. [105]). 

2.7.2. The Far Infrared Background 

The CORE DIRBE (Hauser et al. 1998) and 
FIRAS (Fixsen et al. 1998) experiments detect the 
extragalactic radiation background at A > 100/i. 
The integral for A > 125/i is m = 14 nW m^^ sr^^ 
(Fixsen et al. 1998). Extending the integration 
to 100/i might reasonably be expected to add 2 
nW m~^ sr^^ to this value. The density parameter 
for the sum is 



ripiR = 0.8 ± 0.2 X 10" 



(107) 



This is entry 7.2 in the inventory. It seems to be 
a believable lower bound on the energy density 
in the far infrared. The radiation measurements 
allow room for a comparable amount of energy at 
30 < A < 100 /< (Finkbeincr, Davis & Schlegel 
2000; see also Hauser & Dwek 2001). However, 
this amount of energy shortward of 100/i wave- 
length would distort the TeV 7-ray spectrum from 



the extragalactic source Mrk501 (Quinn 1996; 
Aharonian et al. 1999), by the absorption due 
to e+e^ pair production (Kneiske, Mannheim & 
Hartmann 2002; Konopclko et al. 2003). Thus it 
appears that equation (107) is close to the total 
in the far infrared. 

The optical and infrared light from quasars is 
about one percent of the starlight and far infrared 
backgrounds (§2.5.4). 

2.7.3. A Check: Nuclear Burning 

We may compare the estimate of the present en- 
ergy density in the optical through the far infrared 
against what would be expected from the energy 
stored in the heavy elements (eq. [99]), and what 
would be expected from the picture for stellar evo- 
lution. We comment on the former here and the 
latter in §2.8.1. 

The sum of equations (106) and (107) (entries 
7.2 and 7.3), corrected for the redshift energy loss 
factor in equation (35), is an estimate of the nu- 
clear energy required to produce the observed ra- 
diation, 



riphoton =4.8 ± 1.4x10" 



(108) 



This number may be compared to the nuclear 
binding energy released by nucleosynthesis. The 
estimate in equation (99), reduced by 7% to take 
account of the energy carried away by neutrinos, 
is 



= 5.3 ±1.1 X 10" 



(109) 



^ ^photon, nuc. 

The difference, 10%, is well within our uncertain- 
ties. It might be relevant to note that the differ- 
ence between equations (108) and (109) would be 
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increased if we adopted the estimate of the optical 
background from source counts. 

2.7.4- The X-ray — "f-Ray Background 

Entry 7.4 in Table 1 is the integral of the radia- 
tion background spectrum compiled by Gruber et 
al. (1999) over the energy range 3 keV to 100 GeV. 
The largest contribution to the integral is at pho- 
ton energy ~ 30 keV, but the convergence at high 
energy is slow because the energy per logarithmic 
interval of frequency varies about as ufli, ^ u^^'^. 
The measurement at 3 keV by De Luca & Molendi 
(2004) is about 1.3 times the Gruber et al. value. 
We have not attempted to adjust the entry for this 
more recent result because we do not know its sig- 
nificance for the spectrum at higher energies. 

2.7.5. The Radio Background 

Longair (1995), following Bridle's (1967) study, 
concludes that the brightness temperature of the 
isotropic radio background at 178 MHz, after cor- 
rection for the thermal component, is T = 27±7 K. 
This corresponds to 



6600(1±0.25)(^)-- Jy 



sr 



(110) 



assuming the canonical power law spectrum (see 
also Peacock 1995). Longair also argues that this 
radiation is dominated by galaxies. 

Equation (110) may be compared to the sum of 
the radio source counts in the 9CR radio survey 
at 15 GHz, in the observed range of flux densities, 
0.01 to 1 Jy (Waldram et al. 2003). This sum gives 
Ui, = 5900 Jy sr~^ when scaled to 1 GHz with the 
power index of 0.8. The radio source counts at 
8.4 GHz by Fomalont et al. (2002), summed over 
10 /xJy to 1 Jy, give Uu == 4700 Jy sr^^ at 1 GHz. 
A comparable result is obtained from the count 
at 40 GHZ (as summarized in Figure 13 of Ben- 
nett et al. 2003b). Haarsma & Partridge (1998) 
argue that the integral over the counts likely con- 
verges at 5 '^ 1 /xJy. We conclude that the mea- 
surements are reasonably concordant with equa- 
tion (110), within -0.2 dex. 

The integral of equation (110) slowly diverges 
at short wavelengths. We adopt, as an operational 
definition, a cutoff at A = 1mm (i^ = 300 GHz), 
and we count the contribution at shorter wave- 
lengths as part of the FIR background in our in- 



ventory. There is a natural cutoff at long wave- 
length, at ly ~3 MHz (Simon 1978). Integrating 
eq. (110) over this wavelength range, we obtain 



n 



radio 



5 X 10 



-11 



(111) 



To understand the relation to other energy en- 
tries, we may attempt an alternative estimate by 
directly summing the contributions of known radio 
galaxies. The luminosity function of radio galaxies 
at zero redshift is now reasonably well known by 
virtue of the correlation of a large NRAO VL A Sky 
Survey at 1.4 GHz (NVSS; Condon et al. 1998) 
with galaxies in optical catalogues (Condon, Cot- 
ton & Brodcrick 2002 [UGC vs. NVSS]; Machalski 
& Godlowski 2000 [LCRS vs. NVSS]; Sadler et al. 
2002 [2dF vs. NVSS]). The luminosity function is 
written as the sum of two components, weak radio 
emitters that represent normal galaxies with star 
forming activity, and strong emitters that mostly 
consist of subsets of giant elliptical galaxies and 
AGNs. The former activity is ascribed to electron 
acceleration in Type H (and Ib/c) supernova rem- 
nants. At close to zero redshift the integral over 
the luminosity function gives luminosity densities 



Cg = 1.5 X lO^^W Hz^i Mpc 
Ca = 3.7 X 10^^ W Hz^i Mpc 



-3 



-3 



(112) 



at 1.4 GHz (Condon et al. 2002). The second com- 
ponent, AGNs, is dominated by the most luminous 
sources, and so is subject to sampling fluctuations. 
Haarsma et al. (2000) show that the evolution of 
the flrst component in equation (112) is fast to 
z ^ 1. The evolution they derive is consistent 
with the model for the evolution of the star for- 
mation rate in equation (31). Condon (1992) gives 
the relation between the radio emissivity and the 
star formation rate. After adjustment for the over- 
all constraint from the star formation rate density 
derived from the local Ha luminosity density, as 
discussed in §2.3.2, and our reference IMF, the re- 
lation reads 



'^ - '■' ^ ^0" (ck) 



-0.8 



-2.6 X 10^ 



GHz 



WHz"\ (113) 



per 1 toq yr~^ of star formation. The first 
component represents synchrotron radiation, pri- 
marily from supcrnovae, and the second rep- 
resents bremsstrahlung from HII regions. The 
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bremsstrahlung contribution becomes more im- 
portant above 25 GHz. The integral over the star 
formation history to 2;/ = 3 (with the redshift 
energy loss) and the frequency range v — 3 MHz 
to 300 GHz gives r^radio = 3.3 x 10"^^, where 
the contributions from synchrotron radiation and 
bremsstrahlung are in the proportion 0.4:0.6 in 
the frequency range that concerns us. 

The evolution of the strong radio sources ap- 
pears not to be very fast. The luminosity func- 
tions of Sadler et al. (2002) and Machalski & 
Godlowski (2000) up to z ^ 0.2 suggest a slow 
evolution compared with that of the weak emit- 
ter component, at least at low redshifts. We take 
the evolution factor derived from high redshift 
strong radio galaxies by Willott et al. (2001), 

(1 + z)2-6^(z) with E{z) = [n^{i + z)3 + n^]^/^, 

and we assume the canonical synchrotron spec- 
trum. The integral over time and frequency yields 
f^radio = 2.0 X 10^^^. Thc sum of the two compo- 
nents is 

f^radio - 5 X 10-". (114) 

The consistency with equation (111) suggests that 
this simple model for the radio sources is a useful 
approximation. 

2. 7. 6. Gravitational Radiation 

We consider two sources of gravitational radia- 
tion, the merging of stellar mass binaries and the 
formation of massive black holes in thc centers of 
galaxies. The rate of production of gravitational 
radiation by stellar mass binaries in the Milky 
Way depends on the very uncertain distribution 
of properties of close binary systems that contain 
white dwarfs, neutron stars, and black holes. The 
largest contribution to the mean gravitational ra- 
diation luminosity of the Milky Way is thought to 
be from the merging of binary neutron stars, as 
indicated in Fig. 12 in Schneider et al. (2001). 

Shibata & Uryu (2002) and Faber et al. (2002) 
find that in the merging of a binary neutron star 
system the energy emitted as gravitational radia- 
tion is 

ecwnn ^ 0.005 (115) 

times the mass of thc binary. A major uncertainty 
is the behavior of the system after contact. Gravi- 
tational radiation emitted during the transient for- 
mation of an oscillating neutron star might double 
the number in equation (115). 



If we adopt the estimate of the rate of merg- 
ing of binary neutron stars in the Milky Way, 
7?, « 80 Myr~^, that is based on the three known 
systems that will merge within the Hubble time 
(Kalogera et al. 2004a, b), then we find that 
the product with the radiation energy released 
(eq. [115]), the effective number density of Milky 
Way galaxies (eq. [54]), and 85 Gyr/2 (eqs. [35] 
and [36]) is the gravitational radiation energy den- 
sity r^NS.gw = 10^^. The merger rate is higly un- 
certain: the spread of central values in the models 
in Table 1 in Kalogera et al. (2004b) is lOt?'^- 
On taking account of the upward uncertainty in 
equation (115), we are led to allocate an error of 
a full order of magnitude up or down in entry 7.5. 

The merging rates for neutron star-black hole 
and black hole-black hole binaries are even more 
uncertain. Thc indication from the analysis of 
Schneider et al. (2001) is that mergers involving 
black holes add about 10% to the gravitational ra- 
diation luminosity. 

In the low frequency regime (^ mHz) the lu- 
minosity is dominated by compact white dwarf 
binaries, as shown by Hils, Bender & Webbink 
(1990). Thc integration of the gravitational lu- 
miosity given in their Table 7 with the number of 
close white dwarf binaries in the Milky Way nor- 
malised to 3 X 10^ (after a suggested reduction by 
a factor of ten from the simple theoretical estimate 
of the abundance of these objects) gives the lumi- 
nosity 3 X 10'^^ erg s^^. The product with the 
effective number density of Milky Way galaxies 
and thc effective time span with redshift correc- 
tion gives r^wD.gw ~ 10^^°. This is at the lower 
end of the suggested range of uncertainty in the 
value of 0ns, gw 

Our estimate of the gravitational radiation re- 
leased by the formation of the massive black holes 
in the centers of galaxies is based on the analy- 
sis by Baker et al. (2001), which indicates that 
the gravitational radiation energy released by the 
merging of equal mass black holes is Cngw — 0.03 
times the final black hole mass M, and that 
the frequency of the radiation is on the order of 
0.1M~^, or about 10~^M9 Hz when the mass is 
measured in units of lO^m©. We take the redshift 
at peak production to be z„ ~ 2.5, about the peak 
of quasar activity. Then the present energy den- 
sity of the gravitational radiation from the forma- 
tion of the massive black holes, with mass density 



26 



fin = 10-''•^ is r!„gw ~ e„gw^«/(l + ^«) - 10-^-5. 
This may be compared to the bound on the grav- 
itational wave energy density from pulsar timing 
measurements, ilg < 10^® at frequency ^ 10^^ Hz 
(Kaspi, Taylor & Ryba 1994; Thorsett & Dewey 
1996; Lommen & Backer 2001). The bound is 
about 30% of our estimate of flngw, but there is 
no inconsistency because the expected peak fre- 
quency is considerably larger than the frequency 
where the pulsar constraint is tight. The spec- 
trum of the gravitational radiation from merging 
massive black holes, and the puzzle of how merg- 
ing black holes get close enough that energy loss 
by gravitational radiation drives merging within 
a Hubble time, are discussed by Hughes et al. 
(2001), Jaffc & Backer (2003) and Wyithe & Loeb 
(2003). 

An entry in the inventory for possible sources of 
low frequency gravitational radiation in the early 
universe, as from inflation or phase transitions, 
seems premature. 

2.8. Products of Stellar Evolution 

It is in principle straightforward to compute the 
integrated outputs of stellar evolution — energy, 
neutrinos, helium, and heavy elements — given 
models for the IMF, the star formation history, 
and stellar evolution. Since the details of the re- 
sults of stellar evolution computations are not eas- 
ily assembled, we use approximate estimates by 
procedures similar to those developed in §2.3.1 and 
2.3.2 for the stellar population and its evolution. 
The results add to the checks of consistency of 
our estimates of the stellar production of helium 
and heavy elements and the resulting total energy 
release, and are used to estimate the inventory en- 
tries for the neutrino cosmic energy density. 

2.8.1. Stellar Evolution 

Most stars with masses to < 1 niQ arc still 
on the main sequence. We assume that on av- 
erage 5% of the hydrogen in these subsolar stars 
has been consumed, with energy production ef- 
ficiency 0.0071. Most of the stars with masses 
m > 1 juq have already undergone full evolu- 
tion and left compact remnants, while the fraction 
~ (to/to©) ~^-^ is still on the main sequence. We 
assume that in the latter stars on average 5% of 
the hydrogen has been consumed, as for subsolar 



stars. For the evolved stars we do not attempt 
to follow the details of nuclear burning and mass 
loss. Instead, we adopt estimates of the nuclear 
fuel consumed or mass lost in a few discrete stages 
of evolution, in a similar fashion to the approach 
used in §2.3.1 to tally stellar remnants. 

When the amount of hydrogen consumed in 
stellar burning is 10% of the mass of a star it leaves 
the main sequence. In the model in §2.3.1, stars 
with masses in the range 1 < tti < Stuq eventu- 
ally produce white dwarfs that mainly consist of 
a carbon-oxygen core. In standard stellar evolu- 
tion models, hydrogen burning extends outward 
to a shell after core hydrogen exhaustion, and he- 
lium burning similarly continues in a shell after 
core helium exhaustion. That leaves a CO core 
with the mass given by equation (69). The helium 
layer outside the CO core is thin for m < 2.2?Ti0, 
for which helium ignition takes place as a flash, 
but in stars with masses m > 2.2to,q a signif- 
icant amount of helium is produced outside the 
core, transported by convection to the envelope, 
and liberated. A 5to0 star liberates 0.4 to© of he- 
lium and produces a core with mass 0.85 TOq (e.g., 
Kippenhahn, Thomas & Weigert 1965). We model 
the helium production as a function of the initial 
star mass m as 



mile 



0.23(to- 2.2) + 0.63, 



(116) 



in solar mass units. The energy production in hy- 
drogen shell burning is the product of this mass 
with the post-stellar hydrogen mass fraction, 0.71 
(eq. [87]), and the efficiency factor, 0.0071. Helium 
burning in the core produces energy with efficiency 
factor 0.0010. 

For stars with masses to, > 8?n0 we adopt the 
hehum yield from Table 14.6 of Arnett (1996),^^ 
which we parametrise as 



iTT-Hc — 0.69?Ti — 3.9. 



(117) 



This connects to equation (116) at 8.7?7i0. We 
take the CO core mass as a function of the initial 
stellar mass from Arnett (1996): 

TOCO = 0.28TO- 2.20, to>13to0. (118) 

We use an interpolation of eqs. [69] and [118] for 
stellar masses between 8 and 13to-0. The energy 



^^We prefer cq.[117] which is derived from tabulated values, 
rather than Arnett's pararnetrisation, mjjo = 0.43m — 2.46. 
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Table 5 
Stellar energy production'^ 




stage of stellar evolution 0.08 — Im© 


1 - SniQ 


8 - lOOm© 


sum 


main sequence 0.11 
H shell burning 
core evolution 


0.20 
0.18 
0.05 


0.12 
0.29 
0.04 


0.43 
0.48 
0.09 


sum 0.11 


0.43 


0.46 


1.00 



^Normalised to n = 5.3 X 10 



release is 0.0071 x 0.71 per unit mass for He pro- 
duction, and 0.0014 for CO core formation and the 
further heavy element production. 

The energy output obtained by integration over 
the IMF and PDMF is 



a 



Stella 



E = 5.3 X 10" 



(119) 



The partition into each phase of stellar evolu- 
tion and stellar mass range is given in Table 5, 
where the numbers are normalised to equation 
(119). About 60% of the energy is produced in 
the evolved stages. 

The estimate of the total energy generation in 
equation (119) is in satisfactory agreement with 
our estimate of the nuclear binding energy, JIbe = 
5.7 X 10^^ (eq. [99]), and the energy production 
required to produce our estimate of the present 
radiation energy density, fl-y+i, = (5.1 ± 1.5) x 
10^^ (eq. [108] corrected for neutrino emission, as 
discussed in §2.8.2). That is, our models for stellar 
formation and evolution are consistent with our 
estimates of the accumulation of heavy elements 
and the stellar background radiation. 

In our models for stellar formation and evolu- 
tion the mass of helium that has been liberated to 
interstellar space is 



Arjy w 2.1 X 10" 



(120) 



This may be compared to our estimate from the 
production of heavy elements and the associated 
production of helium, Afly ~ 1-7 ± 1 x 10^** 
(eq. [96]). The heavy element production is cal- 
culated in a similar way. The current wisdom is 
that a neutron star or black hole is left at the end 
of the evolution of a star with mass ?Ti > 8 ttiq , 
and the rest of the mass is returned to the in- 
terstellar medium. On subtracting the remnant 



mass (indicated in Table 2) from the total heavy 
element mass produced, we find that in our model 
the heavy element production is 



O; 



5 X 10" 



(121) 



The more direct estimate in equation (94) is Q'^ ~ 
8 X 10-^ 

We suspect the checks on helium and heavy ele- 
ment production in equations (120) and (121) are 
as successful as could be expected. In particular, 
equation (118) for the CO core mass is not tightly 
constrained. One must consider also the possibil- 
ity that some supernovae do not produce compact 
remnants, as is suggested by the cases of Cas A 
and SN1987A. The maximum heavy element pro- 
duction when there are no remnant neutron stars 
or black holes is three times the value in equa- 
tion (121), and larger than fl'^ (eq. [94]). 

2.8.2. Neutrinos from Stellar Evolution 

We need the fraction /^ of the energy released 
as neutrinos by the various processes of nuclear 
burning. In stars with masses m < 1.4 to©, en- 
ergy generation in hydrogen burning is dominated 



by the slow reaction p + p ^ d + e~^ 



This 



produces neutrinos with mean energy 0.265 MeV, 
which amounts to f^ = 0.020 times the energy 
generated in helium synthesis. In a Solar mass 
star electron capture of ^Be adds 0.005 to the frac- 
tion fi,. The neutrino energy emission fraction is 
larger in higher mass (m > 1.4 to.©) main-sequence 
stars in which the CNO cycle dominates. In this 
process, the neutrinos produced in the Z?"*" decays 
of ^^N and ^^O carry away the energy fraction 
fv = 0.064. The fraction increases to 0.075 for 
TO > 2to0, when the subchain i^N-i'^O-^'^F-^O- 
"'^^N starts dominating and neutrinos are produced 
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by the beta decay of ^^F. This sidcchain also dom- 
inates during shell burning. 

The integral of these factors over the stellar 
IMF normalised to the present-day mass density 
^star = 0.0027 (eq. [27]), together with our pre- 
scription for energy generation in Table 5, yields 
the neutrino energy production, 



f7^^m, ^ 0.31 X 10" 



(122) 



The temperatures and densities that are reached 
up to carbon burning are low enough that there is 
negligible neutrino energy loss from neutrino pair 
production processes. 

The temperatures after carbon burning are high 
enough that the neutrino energy loss dominates, 
that is, /,y ~ 1 (Weaver et al. 1978). Thus we 
may take it that the extra binding energy A(BE) 
of the elements heavier than ^"Ne with respect 
to the binding energy of carbon represents the 
neutrino energy emitted in the late stages of stel- 
lar evolution. On multiplying the heavy element 
mass (the sum of entries 1, 2 and 4 to 9 in Table 
3) by X;(^»/^)A(BE), (where (A(BE)) = 0.0004 
for the solar mix of element abundances, and 
X; Z^jZ = 0.35 is from Grevesse & Sauval (2000), 
we obtain 



17. 



pairs 



= 0.03 X 10 



-6 



(123) 



The sum of equations (122) and (123) is 7% 
of the total energy production (eq. [119]). The 
present energy density of the stellar neutrinos, in 
entry 8.1 in Table 1, is the product of this sum 
with the redshift loss factor ^ 0.5 (eq. [35]). 

2.8.3. White Dwarfs and Neutron Stars 

Most of the gravitational energy liberated in 
white dwarf formation goes to neutrinos, and in 
the latest stage to X-rays. Since the latter is small 
the contribution to the neutrino energy density 
(entry 8.2 in the inventory) is the product of the 
gravitational binding energy in entry 5.3 with the 
redshift factor. 

More than 99% of the energy released in core 
collapse also is carried away by neutrinos. Thus 
we similarly obtain entry 8.3 by multiplying entry 
5.4 by the redshift factor. 



2.9. Cosmic Rays and Magnetic Fields 

This estimate is based on the rate of production 
of cosmic rays in the Milky Way galaxy. We start 
from the radio structure. In their analysis of the 
Haslam et al. (1982) measurements, Beuermann, 
Kanbach, & Berkhuijsen (1985) write the radio 
luminosity distribution at 408 MHz as the sum 
of thin and thick components. We model each 
distribution, in an approximation to their results, 
as 

£(r)=£:(0)e-(l^l/^°''+'''"). (124) 

In the thick component the scale heights in the 
disk and perpendicular to the disk are 



I = 3.1 kpc. 



^cq 



= 1.4 kpc. 



(125) 



scaled to Solar galactocentric distance 8 kpc (see 
Appendix). The total luminosity in the thin com- 
ponent in this model may be neglected, but the 
thin component is a significant contribution to the 
local luminosity density, £q: the local thick com- 
ponent fraction is £0(thick)/f0 = 0.55. We as- 
sume the energy density in cosmic rays is propor- 
tional to the luminosity density at 408 MHz. A 
more detailed model would take account of the 
spatial variation of the magnetic field strength, 
but that will be left for future studies. We nor- 
malize to the local energy density, 1.8 eV cm""^, 
in relativistic cosmic rays (Webber 1998). In these 
approximations the energy of cosmic rays in the 
galaxy is 10^®'^ erg. We take the cosmic ray mean 
life in the galaxy to be 2 x 10^ yr (Garcia-Munoz, 
Mason & Simpson 1977; Yanasak et al. 2001). 
The ratio is a measure of the cosmic ray luminos- 
ity, 10^^'^ erg s^^. If the acceleration of cosmic 
rays is mainly due to shocks of Type II and Ibc 
supernovae,^^ following the conventional wisdom 
(Ginzburg & Syrovatskii 1964), the cosmic ray lu- 
minosity is proportional to the star formation rate. 
Thus the product of this luminosity with the effec- 
tive time (eq. 36) and the redshift factor (eq. [35]) 
is an estimate of the contribution of the Milky Way 
to the present energy in cosmic rays. The product 
with the effective number density of Milky Way 
galaxies (eq. [54]) is our estimate of the present 



^®No type la supernova remnants arc known to give strong 
radio sources (Weiler et al. 1986). This implies that Type 
la supernovae contribute little to the acceleration of cosmic 
rays. 
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energy density in cosmic rays from normal galax- 
ies, 

^CR = 10-*^^ (126) 

The radio luminosity of the Milky Way offers 
a check of this calculation. The Beucrmann ct al. 
(1985) luminosity of the Milky Way, 5.5 x lO^^ 
W Hz-i at 1/ = 408 GHz, scaled to 1.4 GHz by 
the i/~°'® power law (neglecting the small thermal 
bremsstrahlung component), and multiplied by 
the effective number density of Milky Way galaxies 
(eq. [54]), is ^^ = 9 x IQi* W Hz"! Mpc'^, about 
half the Condon et al. (2002) measurement of the 
mean luminosity density of the galaxies (eq. [112]). 
That is, the evidence is that the Milky Way gives 
a reasonably good measure of the synchrotron lu- 
minosity density of the galaxies, and hence of the 
cosmic ray luminosities of galaxies. This is not a 
very direct check of the assumed universality of the 
cosmic ray lifetime in the source galaxy, of course. 

A similar energy may be present in the mag- 
netic field. If the leakage of cosmic rays approxi- 
mated a fluid flow, magnetic field would leak into 
intergalactic space with the cosmic rays. Applica- 
tion of equipartition to equation (126) would sug- 
gest that the cosmic rms magnetic field strength 
is 

BiGM -^ 3 X 10"^ GauB. (127) 

If the estimate of the local ratio of the cosmic ray 
to magnetic field energy density, 8:1 (the magnetic 
field corresponding to 3 /iGaufi), applied to the in- 
tergalactic medium, it would lower equation (127) 
by a factor of three. 

The product of the integrated production of 
supernovae with the characteristic kinetic energy, 
Eke — 1-6 x 10^^ erg, released in a supernova (Ar- 
nett 1996), is an estimate of the integrated kinetic 
energy production per comoving volume. 



n 



SN KE 



10 



-7.3 



(128) 



If the fraction tcr of this energy were placed in rel- 
ativistic intergalactic particles, the present energy 
density in this component, taking account of the 
redshift factor, would be 



n. 



SN cr 



10- 



-7.6, 



(129) 



This is consistent with the picture that a frac- 
tion Ecr ~ 0.2 of the kinetic energy liberated 
by supernovae has been deposited in cosmic rays 
(eq. [126]), which eventually become intergalactic. 



Cosmic rays might gain energy from shocks pro- 
duced by streaming motions of magnetized warm 
plasma in the vicinity of galaxies (Loeb & Wax- 
man 2000). However, since the energy available in 
the plasma (eq. [130] below) is less than the energy 
liberated by supernovae (eq. [126]), this process is 
not likely to substantially affect the energy in cos- 
mic rays. 

What is the AGN contribution to the cosmic 
ray energy density? Though AGNs are not major 
sources at optical wavelengths, their contribution 
to the radio background is comparable to that of 
the larger number of normal galaxies (eq. [112). 
This is an indication of the importance of AGNs 
for high energy processes. If the cosmic ray energy 
production by AGNs were proportional to the pro- 
duction of synchrotron radio emission, the energy 
of cosmic rays from AGN's (including radio ellip- 
tical galaxies) would be comparable to, or even 
larger than, that from the normal galaxies, and 
our estimate in eq. (126) would have to be dou- 
bled. It is important also that AGNs may acceler- 
ate cosmic rays to higher energies than do normal 
galaxies. 

We use equation (126) in the inventory. The 
lower error flag for the entry in Table 1 is based on 
the evidence that the Milky Way gives a fair mea- 
sure of galaxy radio luminosities. The larger upper 
error flag reflects the possibly significant roles of 
the intergalactic magnetic field and AGNs. 

2.10. Kinetic Energy in the IGM 

Intergalactic baryons have peculiar velocities 
that are driven by the gravitational field of the 
dark matter distribution and by nongravitational 
interactions. The former is part of the primeval 
energy in category 4 in the inventory. The latter is 
expected to be most important near the virialized 
regions of galaxies, where the streaming motion 
has been largely transformed by shocks into ther- 
mal energy at the temperature T ^ 2 x 10® K 
associated with the nominal velocity dispersion 
a = 160 km s^^ around L ^ L^, galaxies (eq. [10]). 
This shocked matter may be responsible for the 
O VI absorption lines in the Local Group (Sem- 
bach et al. 2003; Cen et al. 2001). The product 
of the internal energy per unit mass of plasma at 
this temperature with the mass fraction in entry 
3.1a is 

^KE ^ 10"''°. (130) 
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If an appreciable fraction of the kinetic energy 
produced by supernovae (eq. [128]) were deposited 
as kinetic energy in galactic haloes this energy 
would be dissipated by hydrodynamic processes 
rather than the cosmological redshift, but still 
may make a significant addition to equation (130), 
heating the IGM relative to dark matter. 

Baryonic matter well outside the nominal viri- 
alized regions of galaxies and larger systems — at 
distances greater than about r2oo from clusters of 
galaxies and Vy ~ 200 kpc from L ^ L^, galaxies 
— is observed as Lya absorption systems (Pen- 
ton, Stocke & ShuU 2004). The primeval peculiar 
motions of this matter are perturbed by photoion- 
isation that produces kinetic temperatures on the 
order of 10* K. The product of this kinetic energy 
per unit mass with the intergalactic mass frac- 
tion (entry 3.1b) is about one percent of the ki- 
netic energy in the warm intergalactic component 
(eq. [130]). 

We use equation (130) for category 10. We 
caution, however, that a comparable amount of 
kinetic energy may be deposited by supernova 
winds. 

2.11. Electrostatic Energy 

The effect of the electrostatic interaction on the 
binding energies of atomic nuclei relative to free 
protons and electrons is taken into account in cat- 
egory 6, and the electrostatic contribution to the 
binding energies of white dwarfs is (in principle) 
part of category 5. The molecular binding energy 
in objects ranging from dust to asteroids that are 
held together by the electrostatic interaction de- 
serves separate mention. 

The molecular binding energy relative to free 
atoms in condensed matter is roughly 1 eV per 
atom. The product of this mass fraction, ^ 10^^", 
with entry 3.12 is 



Q. 



BE. molecular 



-10" 



(131) 



This is the binding energy of condensed matter 
physics outside strongly self-gravitating systems. 
We refrain from entering it in Table 1 because the 
estimate is so small and uncertain, but we offer for 
comparison the binding energy of the electrons in 
atoms. 

The electrostatic binding energy of the elec- 
trons in an O VI atom is 1.6 keV, and the ad- 



dition of the other five electrons in a neutral oxy- 
gen atom increases the binding energy by only 27 
percent. That is, it is a reasonable approxima- 
tion to ignore the states of ionization of the heavy 
elements and the electrostatic binding energy of 
neutral hydrogen and helium. The sum over the 
heavy element masses in entries 1, 2, 4, 5, and 6 in 
Table 3, weighted by the neutral atomic binding 
energies of 15 cosmically conspicuous elements, is 



VL 



BE, atoms 



-10 



-9.7 



(132) 



This is some six orders of magnitude larger than 
the molecular binding energy in equation (131), 
and comparable in value to the smallest entries in 
Table 1. 

3. Discussion 

The convergence of the cosmological tests from 
redundant constraints on the global cosmologi- 
cal parameters offers a good case for the ACDM 
Friedmann-Lemaitre cosmology as a useful ap- 
proximation to the real world. This is a neces- 
sary condition that our compilation of the energy 
inventory is a meaningful exercise. 

We adopt the working assumption that the 
galaxies are useful tracers of mass, in the sense 
that inventory entries derived under this assump- 
tion are good approximations to reality. This as- 
sumption has been widely questioned. However, 
the recent evidence from weak lensing (eqs. [3] 
to [7]), with the older evidence from the galaxy rel- 
ative velocity dispersion at separations ~ 100 kpc 
to 1 Mpc, and the consistent picture for the viri- 
alized parts of field galaxies (eq. [10]), offer what 
seems to be a reliable case for the use of galaxies as 
mass tracers. The approach certainly is not exact, 
and working out more accurate measures of the 
mass distribution on the scale of galaxies remains 
an important and fascinating challenge. 

Our other conventions for the cosmology are 
less controversial but may have to be adjusted. 
For simplicity, we have adopted a fixed distance 
scale (eq. [1]). We have not presented the scaling 
of inventory entries with the distance scale, which 
can be somewhat complicated; this issue is best re- 
visited when the uncertainty in the distance scale 
is better understood. Our adopted value and for- 
mal uncertainty for the matter density parameter, 
^m (eq. [8]), is in the generally accepted range, but 
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there is reason to suspect that it is somewhat over- 
estimated (§2.1). Again, this is a issue to revisit 
when the constraints have improved. Our value 
for the baryon mass density agrees with what is 
derived from the CMBR anisotropy and the mea- 
surements of the primeval deuterium abundance, 
but it is larger than what is indicated by the 
primeval helium abundance. We have suggested 
in §2.2.2 that the problem may be with the tech- 
nical difficulties of the helium measurements, but 
we will be following the advances in this subject 
with interest. It would not be surprising to see 
the discovery of richer physics in the dark sector, 
and with it an increase in the variety of entries in 
category 1, and possibly also some adjustment of 
other parts of the inventory that depend on the 
dark sector physics. One should also bear in mind 
that the alternative pictures of structure forma- 
tion that were under discussion a decade ago, such 
as cosmic strings, could be operating as subdomi- 
nant perturbations to the ACDM model. Further 
progress in testing the model for structure forma- 
tion will inform our ideas on whether the entries 
that depend on the theory of structure formation 
are likely to require adjustment. Each of these 
issues may be pointing to revisions to the inven- 
tory outside our stated errors. On the other hand, 
the successful network of tests of the cosmology 
and the model for structure formation leads us to 
expect that the general framework presented in 
Table 1 is not likely to change. 

We have argued that 6 ± 1% of the baryons are 
in stars and stellar remnants (eqs. [17] and [27] 
with h — 0.7). This small fraction is analyzed 
in considerable detail in the inventory, in entries 
that are supported by a network of tests. The esti- 
mate in §2.7.1 of the mass density in stars uses the 
optical luminosity density and the stellar-mass-to- 
light ratio. The luminosity density is checked by 
measurements of the optical to near infrared in- 
tergalactic radiation energy density, and by the 
measured galaxy counts (§2.7.1). The comparison 
is not tight, but it shows that the radiation energy 
density likely is known to ±0.3 dex. The conver- 
sion from the luminosity density to the stellar mass 
density depends on the model for the stellar initial 
mass function. The model we use is checked by the 
reasonable agreement with the ratio of the white 
dwarf to subsolar main sequence star densities at 
the low mass end, and with the global Type II su- 



pernova rate at the high mass end (2.3.1, 2.3.2). 
The model for the star formation history is based 
on the time history of the Ha luminosity density, 
which is broadly consistent with other measures 
of the evolution of the star formation rate density 
(subject to the considerable uncertainty in extinc- 
tion in the UV). The model is checked by consis- 
tency with the accumulated mass in stars (after 
correction for mass loss; eq. [34]). Yet another 
check involves the accumulated mass density in 
heavy elements. The release of nuclear binding 
energy in the heavy elements, corrected for the 
loss of radiation energy by the cosmological red- 
shift (eq. [35]), is in satisfactory agreement with 
the present radiation energy density at near op- 
tical and far infrared wavelengths (§2.7.3). This 
test would fail if there were a substantial amount 
of radiation energy at wavelengths 1^ < A < lOO^u; 
improved measurements will be of considerable in- 
terest. This test also depends on the stellar pro- 
duction of helium, which is discussed in §2.2.2, in 
connection with the constraint on the baryon mass 
density from the primeval helium abundance, and 
in §2.8.1, from an analysis of the products of stellar 
evolution. The results seem reasonably consistent. 
The network of checks is complicated, but that is 
what lends credence to the results. 

It is well to pause to consider Arp's (1965) cau- 
tionary remark, that there may be extragalactic 
objects with sizes too small to be readily distin- 
guished from stars — an example is the quasi- 
stellar objects Arp mentions in a note added in 
proof — or with surface brightnesses too low to 
be readily seen against the foreground — an ex- 
ample is the intraclustcr light in the Virgo clus- 
ter (Arp & Bcrtola 1971). The quasar remnants 
could have contained a significant baryon mass, if 
the mass conversion efficiency e„ had been close 
to unity (eq. [73]), but we now know that that 
is difficult to reconcile with the integrated quasar 
energy emission (eq. [83]). There are low surface 
brightness objects (e.g McGaugh, Schombert, & 
Bothun 1995), but the surface brightness of the 
extragalactic sky shows that they cannot largely 
affect our estimate of the mean luminosity density 
(§2.7.1). Substellar objects might be counted as 
part of Arp's cautionary remark, but they are de- 
tectable in nearby galaxies by weak gravitational 
lensing (or MACHOs, as discussed by Alcock et 
al. 2000 and references therein). That is, a half 
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century ago it was not clear that it is feasible to 
establish a fair observational census of the stars. 
Now it appears that the observational conditions 
allow it: the closure of our inventory suggests that 
missing or unknown components cannot be ener- 
getically very significant. 

The baryon number density — excluding 
baryons that may have been sequestered prior to 
light element production — seems to be reliably 
constrained, but the states of most of the baryons 
are not yet observationally well documented. The 
picture that galaxies trace mass leads us to ex- 
pect that about half the dark matter is gathered 
near and within the virialized parts of the galax- 
ies (eq. [13]). Entries 3.1 and 3.2 in Table 1 are 
based on the assumption that the baryons are sim- 
ilarly placed, in the diffuse states observed as hot 
plasma in clusters of galaxies and in warm and hot 
plasma in and around groups of galaxies (eq. [56] ) . 
This has not yet been convincingly observation- 
ally demonstrated. We have not made much use 
of the predicted states of the baryons from numer- 
ical simulations, because we do not know how to 
judge the reliability of the predictions on relatively 
small scales. As the computations and their obser- 
vational checks improve this will become clearer, 
and the results may be expected to improve this 
part of the inventory. 

The advances in observations of the massive 
compact objects in the centers of galaxies allow an 
improved test of the idea that these objects are the 
black hole remnants of quasars and AGNs that are 
powered by gravitational accretion. Our estimates 
of the mass density in black holes and the luminos- 
ity density in quasars require that the quasar mass 
conversion efficiency is e„ ~ 0.02 (eq. [86]), close to 
the efficiency of production of gravitational radi- 
ation (§2.7.6). If the estimate of the quasar lumi- 
nosity density has missed a significant population 
of optically faint AGNs then our value of e„ is bi- 
ased low. The quasar luminosity density is quite 
uncertain, as is the actual value of e„, but it is en- 
couraging that the present result is as close as it 
is to conventional ideas about the quasar emission 
mechanism. 

The eflficiency of conversion of mass into elec- 
tromagnetic radiation in the formation of stellar 
mass black holes is constrained to be less than 
a few percent by the observed energy density in 
electromagnetic radiation and neutrinos. The in- 



terpretation of this constraint awaits development 
of the theory of the astrophysical processes of for- 
mation of stellar mass black holes. Our estimate of 
the gravitational binding energy released by core 
collapse supernovae (entry 8.3) is comparable to 
the energy released in the formation of the mas- 
sive black hole quasar remnants, and it is not very 
much less than the upper bound from neutrino 
detectors. Progress here will be followed with in- 
terest. 

The largest energy densities in the inventory 
are in dark matter and dark energy. Any em- 
pirical hints to the natures of these components 
will be followed with great interest. One possibil- 
ity is that the observations of microlensing toward 
the Magellanic Clouds have detected massive com- 
pact dark objects (Alcock et al. 2000; Afonso et 
al. 2003). We have adopted the conservative in- 
terpretation, that the lenses are baryonic — stars 
and stellar remnants — with standard mass func- 
tions, and that the high lensing rate is an accident 
of the distributions of source and lensing stars. 
A demonstration of the apparently simpler inter- 
pretation, that the lenses contribute some 20% of 
the mass of the halo of the Milky Way, and are 
baryonic, perhaps white dwarfs, would consider- 
ably upset our considerations of the stellar baryon 
budget; a demonstration that the lenses are non- 
baryonic would profoundly affect thinking about 
dark matter and the early universe. 

Though there are large gaps in our understand- 
ing of the cosmic energy inventory, most notably 
in the dark sector, there also is an impressively 
broad observational basis for many parts of the 
inventory. Ten years ago the basis was consider- 
ably smaller, and we can be sure that ten years 
from now the observations will allow more entries 
and a richer network of crosschecks of the inter- 
pretations. Progress will be uneven, of course, as 
one sees by considering the prospects for further 
checks of the neutrino energy inventory. With the 
recent demonstration that the solar neutrino lumi- 
nosity agrees with standard physics augmented by 
neutrino masses and mixing angles (Bahcall 2003 
and references herein) we may now be confident 
that the theory of neutrino production by stel- 
lar nuclear burning is reliably tested even in the 
absence of a direct observation of stellar neutri- 
nos from remote stars. The neutrinos produced 
by plasma processes in white dwarf formation are. 
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while an enormous energy density, too soft for de- 
tection by any conceivable method: estimates of 
this component of the inventory will have to con- 
tinue to rely on the theory. Neutrinos released in 
stellar core collapse have been detected (Hirata et 
al. 1987; Bionta et al. 1987), and the integrated 
background of relic neutrinos from this process is 
expected to be within experimental reach. The 
detection would complete another valuable check 
of the cosmic energy transactions. 

Parts of our inventory can be improved on the 
basis of what is available within existing compu- 
tations, or could be readily developed from exist- 
ing work, but that we could not readily assemble. 
In particular, the theoretical analyses of stellar 
remnants and the other products of stellar evo- 
lution (§§2.3.1 and 2.8), with special attention to 
the white dwarfs that store so much of the heavy 
elements, will have to be done better. Another 
example is the heavy element abundances, which 
can be classified by element type, as oxygen — 
which has a specific spectroscopic importance — 
and the a-elements and the iron group, and by en- 
vironment, as white dwarfs, main sequence stars, 
and interstellar and intergalactic space. 

It is also worth stating explicitly that we may 
expect that, as in the past, there will be quali- 
tatively new experimental or observational devel- 
opments — perhaps the detection of gravitational 
radiation, or the identification of dark matter par- 
ticles — that substantially affect the directions of 
research and the development of the cosmic energy 
inventory. 



demonstrate that a considerable part of the in- 
ventory is a believable approximation. Continued 
advances in the observational and theoretical basis 
for the inventory surely will yield unexpected re- 
visions and additions; we are attempting to draw 
large conclusions from limited observations of an 
exceedingly complex universe. However, the big 
surprise at the moment is that it is now possible 
to find an inventory with observational support for 
the largest ^ 40 forms of energy. 
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4. Concluding Remarks 

A half century ago we had a largely conjectural 
picture of the large-scale structure of the physi- 
cal universe, of its material contents, and of the 
main process that drive transformations among 
the states of matter and radiation. Now our world 
picture has a substantial basis in experiments and 
observations and the attendant theories. This has 
been an evolutionary process: none the entries 
in the inventory in Table 1 requires a substantial 
departure from ideas that are under discussion, 
and in many cases have been in the literature for 
decades. The important new development is that 
we now have observational support for the many 
entries in Table 1, and a network of tests that 
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A. A Few Astrophysical Quantities for the Milky Way 

The current best value for the distance to the Galactic center is 

i?o = 7.94 ± 0.42 kpc, (Al) 

from measurements of the orbit — proper motion and redshift — of a star close to Sgr A* (Eisenhauer et 
al. 2003). The product of this distance with the measurement of the proper motion of Sgr A* by Reid & 
Brunthaler (2004) is 241 km s^^. This velocity corrected for the solar proper motion is our adopted value 
of the circular speed of the Milky Way at the Solar circle, 

Oo = 234±13 kms"\ (A2) 

Traditionally, one invokes the Oort model to infer the circular speed, but Oiling & Dehnen (2003) find that 
the current best result depends substantially on the stellar populations used in the analysis. 

The circular velocity yields an estimate of the optical luminosity of the Milky Way via the TuUy-Fisher 
relation. To convert the optically measured maximum rotation velocity Vrot to the 20% line width of HI 
measurements 14^20 7 we use the relation 

W^2o/2 = (0.93 ± O.G2)K.ot + (27 ± 1.5)km s~\ (A3) 

derived from the data given by Mathewson & Ford (1996). From the Galactic rotation curve from HI and CO 
terminal velocities (using Fig. 9.17 of Binney & Merrifield 1998), we estimate that the maximum rotation 
velocity of the Milky Way is Knax = 1.036o = 241 ± 13 km s^^. Identifying Vmax — K-ot, we find 

14^20 (MW) == 502 ± 30 km s"\ (A4) 

From the COBE-DIRBE integrated light for the Milky Way in the near infrared (J, K, L) wavelength 
bands, Malhotra et al. (1996) showed that the Milky Way is on the standard TuUy-Fisher relation in these 
wave length bands, slightly on the brighter side but within the dispersion around the relation. Assuming 
that the Milky Way has the standard color of spiral galaxies, we may expect that the luminosity in the B 
band also satisfies the TuUy-Fisher relation for the B band. Using the Sakai et al. (2000) local calibration, 

Mb^ = -19.80 - 7.97(log 14^20 - 2.5), (A5) 

we find 

Lb= (MW) = 1O1"-^4±o-19L0, (A6) 

where we use Bq = 5.46, and the error is dominated by the intrinsic dispersion (0.43 mag) of the B-band 
TuUy-Fisher relation. The superscript c stands for extinction-free quantities. 

An alternative estimate of the B-band luminosity uses the colour transformation applied to the DIRBE 
integrated K-band luminosity of the Milky Way. A comparison of the samples of Malhotra et al.(1996) and 
Sakai et al. (2000) indicates B^-K"^ = 2.84 ±0.10 for typical Sbc galaxies. Thus /i:'=(MW) = -23.95 ±0.25 
yields B"^ = -21.1 ±0.3, or 

Lbc(MW) = 1O1"-^3"-13L0. (A7) 

Combining the two, our recommended value is 

Lbo(MW) =4.6±1.4x 10^"Lq. (A8) 

A value often quoted in the literature, i'^ = 2.3 ± 0.6 x 10^°Lq (van den Bergh 1988), is 1.5 standard 
deviations below equation (A8). 

The luminosity obtained here is the extinction-free value. When it is compared with L* that appears in 
the luminosity function, we must take account of internal extinction. For consistency we use the average 
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value of the internal extinction correction in the Sakai et al. (2000) local calibration, Ag*^ = 0.44 mag, or 
an attenuation factor 1.5. This means that Lb(MW) — I.4OL3 at h = 0.7, and that the effective number 
density of Milky Way galaxies is the value in equation (54), jCb/Lb{MW) = O.OlSh^ Mpc"^. 



36 



REFERENCES 

Abazajian, K., et al. 2003, AJ, 126, 2081 
Afonso, C. ct al. 2003, A&A, 400, 951 
Aharonian, F. A., et al. 1999, A&A, 349, 11 

Alcock, C. et al. (The MACHO Collaboration) 

2000, ApJ, 542, 281 

AUende Prieto, C, Lambert, D. L., & Asplund, 
M. 2001, ApJ, 556, L63 

Arnett, D. 1996, Supernovae and Nucleosynthe- 
sis: An Investigation of the History of Matter, 
from the Big Bang to the Present (Princeton: 
Princeton University Press) 

Arp, H. 1965, ApJ, 142, 402 

Arp, H. & Bertola, F. 1971, ApJ, 163, 195 

Bahcah, J. N. 2003, astro-ph/G31G03G 

Bahcall, J. N. & Peha-Garay, C. 2003, JHEP 0311, 
004 

Bahcall, J. N. & Pinsonncauh, M. H. 2004, Phys. 
Rev. Lett. 92, 12130 

Bahcall, J. N., Pinsonncault, M. H., & Basu, S. 

2001, ApJ, 555, 990 (BP2000) 

Bahcall, N. A. & Bode, P. 2003, ApJ, 588, LI 

Bahcall, N. A. & Cen, R. 1993, ApJ, 407, L49 

Bahcall, N. A., Lubin, L. M., & Dorman, V. 1995, 
ApJ, 447, L8 

Bahcall, N. A., et al. 2003, ApJ, 585, 182 

Baker, J., Briigmann, B., Campanelli, M. Lousto, 
C. O. & Takahashi, R. 2001, Phys. Rev. Lett. 
87, 121103 

Basu, B. & Lynden-Bell, D. 1990, QJRAS, 31, 359 

Bennett, C. L. et al. 2003a, ApJS, 148, 1 

Bennett, C. L., et al. 2003b, ApJS, 148, 97 

Bernstein, R. A., Freedman, W. L., & Madore, 
B. F. 2002, ApJ, 571, 56 

Bcuermann, K., Kanbach, C, & Berkhuijscn, 
E. M. 1985, A&A, 153, 17 



Binney , J. & Merrifield, M. 1998, Galactic 
Astronomy (Princeton: Princeton University 
Press) 

Bionta, R. M., Blcwitt, C, Bratton, C. B., 
Caspere, D., & Ciocio, A. 1987, Physical Re- 
view Letters, 58, 1494 

Blanton, M. R., et al. 2003, ApJ, 592, 819 

Bosma, A. 1981, AJ, 86, 1825 

Branch, D. 1992, ApJ, 392, 35 

Bridle, A. H. 1967, MNRAS, 136, 219 

Bruzual, G. & Chariot, S. 2003, MNRAS, 344, 
1000 

Burgasser, A. J., Kirkpatrick, J. D., McElwain, 
M. W., Cutri, R. M., & Skrutskie, M. F. 2003, 
AJ, 125, 850 

Burgasser, A. J., McElwain, M. W., Kirkpatrick, 
J. D., Cruz, K. L., Tinney, C. G., & Reid, I. N. 
2004, AJ, 127, 2856 

Burrows, A. 1990, in Supernovae, ed. by Petschck, 
A. G. (Berlin: Springer), 143 



Burrows, A., Hubbard, W. B., Lunine, J. I., & 
Licbert, J. 2001, Reviews of Modern Physics, 
73, 719 

Calura, F. & Mattcucci, F. 2004, MNRAS, 350, 
351 

Cappellaro, E., Evans, R., & Turatto, M. 1999, 
A&A, 351, 459 

Cappellaro, E., Turatto, M., Tsvetkov, D. Y., Bar- 
tunov, O. S., PoUas, C, Evans, R., & Hamuy, 
M. 1997, A&A, 322, 431 

Cen, R., Tripp, T. M., Ostriker, J. P., & Jenkins, 
E. B. 2001, ApJ, 559, L5 

Chabrier, G. 2003, PASP, 115, 763 

Chen, H., Lanzetta, K. M., Webb, J. K., & Bar- 
cons, X. 2001, ApJ, 559, 654 

Chiosi, C, BerteUi, G., & Bressan, A. 1992, 
ARA&A, 30, 235 

Chokshi, A. & Turner, E. L. 1992, MNRAS, 259, 
421 



37 



Churchill, C. W., Vogt, S. S., & Charlton, J. C. 
2003, AJ, 125, 98 

Claver, C. F., Liebert, J., Bergeron, P., & Koester, 
D. 2001, ApJ, 563, 987 

Cole, S. et al. 2001, MNRAS, 326, 255 

Colless, M., et al. 2001, MNRAS, 328, 1039 

Condon, J. J. 1992, ARA&A, 30, 575 

Condon, J. J., Cotton, W. D., & Broderick, J. J. 
2002, AJ, 124, 675 

Condon, J. J., Cotton, W. D., Greiscn, E. W., Yin, 
Q. F., Parley, R. A., Taylor, G. B., & Broderick, 
J. J. 1998, AJ, 115, 1693 

Cox, A. N. (ed) 2000, Allen's Astrophysical Quan- 
tities, 4th ed. (Berlin: Springer) 

Cox, J. P. & GiuUi, R. T. 1968, Principles of Stel- 
lar Structure, Vol. 2 (New York: Gordon and 
Breach) 

Groom, S. M., Smith, R. J., Boyle, B. J., Shanks, 
T., Miller, L & Outram 2004, astro-ph/0403040 

Davis, M. & Peebles, P. J. E. 1983, ApJ, 267, 465 

De Luca, A. & Molendi, S. 2004, A&A, 419, 837 

de Sitter, W. 1917, MNRAS, 78, 3 

de Vaucouleurs, G. & Olson, D. W. 1982, ApJ, 
256, 346 

Draine, B. T. 2003, ARA&A, 41, 241 

Eguchi, K. et al. 2003, Phys. Rev. Lett. 90, 021802 

Einstein, A. 1917, Sitzungsber Preufi Akad Wiss, 
p. 142 [English translation in The Principle of 
Relativity (Dover, New York, 1952), p. 177] 

Eisenhauer, F., Schodel, R., Genzel, R., Ott, T., 
Tecza, M., Abutcr, R., Eckart, A., & Alexan- 
der, T. 2003, ApJ, 597, L121 

Elvis, M., et al. 1994, ApJS, 95, 1 

Ezer, D. & Cameron, A. G. W. 1967, Canadian 
Journal of Physics, 45, 3429 

Faber, J. A., Grandclcmcnt, P., Rasio, F. A., & 
Taniguchi, K. 2002, Physical Review Letters, 
89, 231102 



Fan, X., et al. 2001, AJ, 121, 54 

Fang, T., Marshall, H. L., Lee, J. C., Davis, D. S. 
& Canizares, C. R. 2002, ApJ, 572, L127 

Fang, T., Sembach, K. R., & Canizares, C. R. 
2003, ApJ, 586, L49 

Ferrarese, L. 2002, Proceedings, Second KIAS As- 
trophysics Workshop, Seoul, 2001, September 
3-7, Singapore: World Scientific, edited by Lee, 
C.-H., 3 

Finkbeiner, D. P., Davis, M., & Schlegcl, D. J. 
2000, ApJ, 544, 81 

Fixsen, D. J., Cheng, E. S., Gales, J. M., Mather, 
J. C., Shafcr, R. A., & Wright, E. L. 1996, ApJ, 
473, 576 

Fixsen, D. J., Dwek, E., Mather, J. C, Bennett, 
C. L., & Shafer, R. A. 1998, ApJ, 508, 123 

Fomalont, E. B., Kellermann, K. I., Partridge, 
R. B., Windhorst, R. A., & Richards, E. A. 
2002, AJ, 123, 2402 

Freedman, W. L., et al. 2001, ApJ, 553, 47 

Fukazawa, Y., Makishima, K., Tamura, T., Ezawa, 
H., Xu, H., Ikebe, Y., Kikuchi, K., & Ohashi, 
T. 1998, PASJ, 50, 187 

Fukuda, Y., et al. 1998, Physical Review Letters, 
81, 1562 

Fukugita, M. 2001, NATO ASIC Proc. 565: Struc- 
ture Formation in the Universe, 93 

Fukugita, M., Hogan, C. J., & Peebles, P. J. E. 
1998, ApJ, 503, 518 (FHP) 

Fukugita, M. & Kawasaki, M. 2003, MNRAS, 340, 
L7 

Gallego, J., Zaniorano, J., Aragon-Salamanca, A., 
& Rego, M. 1995, ApJ, 455, LI 

Gal- Yam, A. & Maoz, D. 2004, MNRAS, 347, 942 

Gao, L., Loeb, A., Peebles, P. J. E., White, S. D. 
M. & Jenkins, A. 2003, astro-ph/03 12499 

Garcia-Munoz, M., Mason, G. M., & Simpson, 
J. A. 1977, ApJ, 217, 859 



38 



Ginzburg, V. I. & Syrovatskii, S. I. 1964, The Ori- 
gin of Cosmic Rays (Oxford: Pcrgamon) 

Glazebrook, K., Blake, C, Economou, F., Lilly, 
S., & CoUess, M. 1999, MNRAS, 306, 843 

Glazebrook, K., et al. 2003, ApJ, 587, 55 

Gott, J. R., Schramm, D. N., Tinsley, B. M., & 
Gunn, J. E. 1974, ApJ, 194, 543 

Gebhardt, K., et al. 2000, ApJ, 539, L13 

Gould, A., BahcaU, J. N., & Flynn, C. 1996, ApJ, 
465, 759 

Grevesse, N. & Sauval, A. J. 2000, Origin of Ele- 
ments in the Solar System, Implications of Post- 
1957 Observations, Proceedings of the Interna- 
tional Symposium, ed. O. Manuel. (Dordrecht: 
Kluwer Academic Publishers), 261 

Gruber, D. E., Matteson, J. L., Peterson, L. E., & 
Jung, G. V. 1999, ApJ, 520, 124 

Gush, H. P., Halpern, M., & Wishnow, E. H. 1990, 
Physical Review Letters, 65, 537 

Haarsma, D. B. & Partridge, R. B. 1998, ApJ, 503, 
L5 

Haarsma, D. B., Partridge, R. B., Windhorst, 
R. A., & Richards, E. A. 2000, ApJ, 544, 641 

Harris, H. C., et al. 2003, AJ, 126, 1023 

Haslam, C. G. T., Stoffel, H., Salter, C. J., & Wil- 
son, W. E. 1982, A&AS, 47, 1 

Hauser, M. G. & Dwek, E. 2001, ARA&A, 39, 249 

Hauser, M. G., et al. 1998, ApJ, 508, 25 

Hawkins, E. et al. 2003, MNRAS 346, 78 

Hayashi, C., & Nakano, T. 1963, Prog. Theor. 
Phys. 30, 460 

Heavens, A., Panter, B., Jimenez, R., & Dunlop, 
J. 2004, Nature, 428, 625 

Heger, A., Woosley, S. E., Fryer, C. L., & Langer, 
N. 2003, From Twilight to Highlight: The 
Physics of Supernovae, 3 (astro-ph/0211062) 

Henry, R. C. 1999, ApJ, 516, L49 



Hils, D., Bender, P. L., & Webbink, R. F. 1990, 
ApJ, 360, 75 

Hirata, K., Kajita, T., Koshiba, M., Nakahata, 
M., & Oyama, Y. 1987, Physical Review Let- 
ters, 58, 1490 

Hubble, E. P. 1929, Proc Nat Acad Sci, 15, 168 

Hubble, E. P. 1934, ApJ, 79, 8 

Hubble, E. P. 1936, The Realm of the Nebulae 
(New Haven: Yale University Press) 

Huchra, J., Jarrett, T., Skrustkie, M., Cutri, R., 
Schneider, S., Macri, L., Steining, R., & Mader, 
J. 2003, lAU Symposium, 216, 174 

Hughes, S. A., Marka, S., Bender, P. L. & Hogan, 
C. J. 2001, in Proceedings of Snowmass 2001, 
APS/DPF/DPB Summer Study on the Future 
of Particle Physics, cConf C010630: p402 (also 
astro-ph/01 10349) 

Irvine, W. M. 1961, Doctoral Dissertation, Har- 
vard University 

Ivezic, Z., et al. 2000, AJ, 120, 963 

Izotov, Y. I. & Thuan, T. X. 2004, ApJ, 602, 200 

Jaffe, A. H. & Backer, D. C. 2003, ApJ, 583, 616 

Janka, H.-T. & HiUebrandt, W. 1989, A&A, 224, 
49 

Jimenez, R., Flynn, C., MacDonald, J., & Gibson, 
B. K. 2003, Science, 299, 1552 

Kalogera, V., et al. 2004a, ApJ, 601, L179 

Kalogera, V., et al. 2004b, astro-ph/0312101, v. 3 

Kameda, J. et al. 2001, in Proc. International 
Cosmic Ray Conference, Hamburg 2001, ed. 
K.-H. Kampert, G. Heizelmann & C. Spiering 
(Katlenburg-Lindau: Copernicus Gesellschaft, 
2001) Vol. 3, 1057 

Kaspi, V. M., Taylor, J. H., & Ryba, M. F. 1994, 
ApJ, 428, 713 

Kauffmann, G. et al. 2003, MNRAS, 341, 33 

Kennicutt, R. C. 1983, ApJ, 272, 54 

Keres, D., Yun, M. S., & Young, J. S. 2003, ApJ, 
582, 659 



39 



Kippenhahn, R., Thomas, H. C, & Weigert, A. 
1965, Zcits. ffir Astrophysik, 61, 241 

Kirkman, D., Tytlcr, D., Suzuki, N., O'Meara, 
J. M., & Lubin, D. 2003, ApJS, 149, 1 

Kneiske, T. M., Mannheim, K., & Hartmann, 
D. H. 2002, A&A, 386, 1 

Kobulnicky, H. A. & Zaritsky, D. 1999, ApJ, 511, 
118 

Kochanek, C. S. et al. 2001, ApJ, 560, 566 

Konopelko, A., Mastichiadis, A., Kirk, J., de 
Jager, O. C, & Stacker, F. W. 2003, ApJ, 597, 
851 

Kormendy, J. & Richstone, D. 1995, ARA&A, 33, 
581 

Kroupa, P. 2001, MNRAS, 322, 231 

Kumar, S. S. 1963, ApJ, 137, 1121 

LaRoque, S. J., et al. 2003, ApJ, 583, 559 

Layzer, D. 1963, ApJ, 138, 174 

Lilly, S. J., Le Fevre, O., Hammer, F., & Cramp- 
ton, D. 1996, ApJ, 460, LI 

Loeb, A. & Waxman, E. 2000, Nature, 405, 156 

Lommen, A. N. & Backer, D. C. 2001, BAAS, 33, 
1347 

Longair, M. 1995, Proceedings of the Extragalac- 
tic Background Radiation meeting, Baltimore, 
1993 May 18-20, Space Telescope Science Insti- 
tute Symposium Series, meeting, ed. Calzetti, 
D., Livio, M., & Madau, P. (Cambridge: Cam- 
bridge University Press), 223 

Machalski, J. & Godlowski, W. 2000, A&A, 360, 
463 

Madau, P., della Valle, M., & Panagia, N. 1998, 
MNRAS, 297, L17 

Madau, P., Ferguson, H. C, Dickinson, M. E., 
Giavahsco, M., Steidel, C. C, & Fruchter, A. 
1996, MNRAS, 283, 1388 

Madau, P. & Pozzetti, L. 2000, MNRAS, 312, L9 

Malhotra, S., Spergel, D. N., Rhoads, J. E., & Li, 
J. 1996, ApJ, 473, 687 



Marcy, G. W. & Butler, R. P. 2000, PASP, 112, 
137 

Mather, J. C, et al. 1990, ApJ, 354, L37 



Mather, J. C, Fixsen, D. J., Shafer, R. A., Mosier, 
C, & Wilkinson, D. T. 1999, ApJ, 512, 511 

Mathewson, D. S. & Ford, V. L. 1996, ApJS, 107, 
97 

Mattila, K. 1990, lAU Symp. 139: The Galactic 
and Extragalactic Background Radiation, 139, 

257 

McGaugh, S. S., Schombert, J. M., & Bothun, 
G. D. 1995, AJ, 109, 2019 

McKay, T. A. et al. 2001, astro-ph/0108013 

Merritt, D. & Ferrarese, L. 2001, ApJ, 547, 140 

Molster, F. J., Waters, L. B. F. M., Tielens, 
A. G. G. M., Koike, C, & Chihara, H. 2002, 
A&A, 382, 241 

Mushotzky, R. F., Cowie, L. L., Barger, A. J., & 
Arnaud, K. A. 2000, Nature, 404, 459 

Mushotzky, R. F. & Loewenstein, M. 1997, ApJ, 
481, L63 

Nakamura, O., Fukugita, M., Brinkmann, J., & 
Schneider, D. P. 2004, AJ, 127, 2511 

Nakamura, O., Fukugita, M., Yasuda, N., Love- 
day, J., Brinkmann, J., Schneider, D. P., Shi- 
masaku, K., & SubbaRao, M. 2003, AJ, 125, 
1682 

Nice, D. J., Splaver, E. M., & Stairs, I. H. 2003, 
ASP Conf. Ser. 302: Radio Pulsars, 75 

Nomoto, K. 1984, ApJ, 277, 791 

Oiling, R. P. & Dehnen, W. 2003, ApJ, 599, 275 

Pagel, B. E. J. & Portinari, L. 1998, MNRAS, 298, 
747 

Pain, R., et al. 2002, ApJ, 577, 120 

Partridge, R. B. & Peebles, P. J. E. 1967, ApJ, 
148, 377 



40 



Peacock, J. A. 1995, Proceedings of the Extra- 
galactic Background Radiation meeting, Bal- 
timore, 1993 May 18-20, Space Telescope Sci- 
ence Institute Symposium Series, meeting, ed. 
by Calzetti, D., Livio, M., & Madau, P. (Cam- 
bridge: Cambridge University Press), 237 

Peebles, P. J. E. 1986, Nature, 321, 27 

Peebles, P. J. E. & Partridge, R. B. 1967, ApJ, 
148, 713 

Peebles, P. J. & Ratra, B. 2003, Reviews of Mod- 
ern Physics, 75, 559 

Peimbert, M., Peimbert, A., & Ruiz, M. T. 2000, 
ApJ, 541, 688 

Penton, S. V., Stocke, J. T. & ShuU, J. M. 2004, 
astro-ph/0401036 

Quinn, J., et al. 1996, ApJ, 456, L83 

Reid, I. N., Gizis, J. E., & Hawlcy, S. L. 2002, AJ, 
124, 2721 

Reid, M. J. & Brunthaler, A. 2004, astro- 
ph/0408107 

Reimers, D. & Koester, D. 1982, A&A, 116, 341 

Reiprich, T. H. & Bohringer, H. 2002, ApJ, 567, 
716 

Richter, P., Savage, B. D., Sembach, K. R., Tripp, 
T. M., & Jenkins, E. B. 2003, ASSL Vol. 281: 
The IGM/Galaxy Connection. The Distribu- 
tion of Baryons at z=0, 85 

Riess, A. et al. 2004, ApJ, 607, 665 

Rudnick, G., et al. 2003, ApJ, 599, 847 

Sadler, E. M., et al. 2002, MNRAS, 329, 227 

Sakai, S., et al. 2000, ApJ, 529, 698 

Salpeter, E. E. 1955, ApJ, 121, 161 

Salucci, P. & Persic, M. 1999, MNRAS, 309, 923 

Salucci, P., Szuszkiewicz, E., Monaco, P., & 
Danese, L. 1999, MNRAS, 307, 637 

Scalo, J. M. 1986, Fundamentals of Cosmic 
Physics, 11, 1 



Schwarzschild, M. 1958, Structure and Evolution 
of the Stars (Princeton: Princeton University 
Press) 

Searle, L., Sargent, W. L. W., & Bagnuolo, W. G. 
1973, ApJ, 179, 427 

Schneider, R., Ferrari, V., Matarrese, S., & Porte- 
gies Zwart, S. F. 2001, MNRAS, 324, 797 

Sembach, K. R., et al. 2003, ApJS, 146, 165 

Shapiro, S. L. & Teukolsky, S. A. 1983, Black 
Holes, White Dwarfs, and Neutron Stars: The 
Physics of Compact Objects (New York: Wiley- 
Interscience) 

Sheth, R. K., et al. 2003, ApJ, 594, 225 

Shibata, M. & Uryu, K. 2002, Progress of Theo- 
retical Physics, 107, 265 

Shull, J. M. 2003, ASSL Vol. 281: The 
IGM/Galaxy Connection. The Distribution of 
Baryons at z=0, 1 

Shull, J. M., Tumlinson, J., & Giroux, M. L. 2003, 
ApJ, 594, L107 

Simon, A. J. B. 1978, MNRAS, 184, 537 

Soltan, A. 1982, MNRAS, 200, 115 

Spergel, D. N. et al. 2003, ApJS, 148, 175 

Steidel, C. C, Adelberger, K. L., Giavalisco, M., 
Dickinson, M., & Pcttini, M. 1999, ApJ, 519, 1 

Steidel, C. C, Dickinson, M., & Pcrsson, S. E. 
1994, ApJ, 437, L75 

Tammann, G. A., Loeffler, W., & Schroeder, A. 
1994, ApJS, 92, 487 

Tegmark, M. et al. 2004a, Physical Review D69, 
103501 

Tegmark, M. et al. 2004b, ApJ, 606, 702 

Thorsett, S. E. & Chakrabarty, D. 1999, ApJ, 512, 

288 

Thorsett, S. E. & Dewey, R, J. 1996, Phys Rev 
D53, 3468 

Tinslcy, B. M. 1973, A&A, 24, 89 

Tremaine, S., et al. 2002, ApJ, 574, 740 



41 



Tripp, T. M., Savage, B. D., & Jenkins, E. B. 2000, 
ApJ, 534, LI 

Trujillo, C. A., Jewitt, D. C, & Luu, J. X. 2001, 
AJ, 122, 457 

van den Bergh, S. 1988, Comments on Astro- 
physics, 12, 131 

van den Bergh, S. & McClurc, R. D. 1994, ApJ, 
425, 205 

Vennes, S., Smith, R. J., Boyle, B. J., Croom, 
S. M., Kawka, A., Shanks^ T., MiUer, L., & 
Loaring, N. 2002, MNRAS, 335, 673 

Verde, L., et aL 2002, MNRAS, 335, 432 

Waldram, E. M., Pooley, G. G., Grainge, K. J. B., 
Jones, M. E., Saunders, R. D. E., Scott, P. F., 
& Taylor, A. C. 2003, MNRAS, 342, 915 

Weaver, T. A., Zimmerman, G. B., & Woosley, 
S. E. 1978, ApJ, 225, 1021 

Webber, W. R. 1998, ApJ, 506, 329 

Weidemann, V. 2000, A&A, 363, 647 

Weiler, K. W., Sramek, R. A., Panagia, N., van 
der Hulst, J. M., & Salvati, M. 1986, ApJ, 301, 
790 

Whelan, J. & Iben, I. J. 1973, ApJ, 186, 1007 

White, D. A. 2000, MNRAS, 312, 663 

Willott, C. J., Rawlings, S., BlundcU, K. M., Lacy, 
M., & Eales, S. A. 2001, MNRAS, 322, 536 

Worsley, M. A. et al. 2004, astro-ph/0404273 

Wright, E. L. & Reese, E. D. 2000, ApJ, 545, 43 

Wyithe, J. S. B. and Locb, A. 2003, ApJ, 590,691 

Yanasak, N. E., et al. 2001, ApJ, 563, 768 

York, D. G., et al. 2000, AJ, 120, 1579 

Zhang, P., Pen. W.-L. & Trac, H. 2004, astro- 
ph/0402115 

Zwaan, M. A. et al. 2003, AJ, 125, 2842 

Zwicky, F. 1933, Hclv. Phys. Acta, 6, 110 



This 2-coluimi preprint was prepared with the AAS lATJiiX 
macros v5.2. 



42 



